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This  is  the  final  technical  report  prepared  by  Calspan 
Corporation  on  a  portion  of  a  multi-phase  program  sponsored  by  the 
Air  Force  Aero -Propuls ion  Laboratory,  Air  Force  Systems  Command, 
Wright -Patterson  AFB,  Ohio  under  Contract  F33615-76-C-2092 . 


The  work  reported  herein  was  performed  as  part  of  Phases  I  and  II 
of  Project  3066,  "Investigation  of  Rotating  Stall  and  Turbine  Heat 
Transfer  in  Axial  Flow  Turbomachinery,"  with  Mr.  Marvin  A.  Stibich, 
AFAPL/TBC ,  as  Project  Engineer.  Dr.  Gary  R.  Ludwig  of  the  Aero¬ 
dynamic  Research  Department  of  Calspan  Corporation,  had  overall 
technical  responsibility  for  the  program,  and  carried  out  the 
experimental  measurements.  The  theoretical  studies  of  rotating  stall 
were  performed  by  Dr.  Joseph  P.  Nenni.  Mr.  John  Nemeth  designed  the 
experimental  apparatus  and  assisted  in  the  experiments . 
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SECTION  I 
INTRODUCTION 


The  optimum  performance  of  a  turbo-propulsion  system  is  usually 
achieved  when  the  compressor  is  operating  near  its  maximum  pressure  ratio. 
However,  this  optimum  is  generally  not  attainable  because  it  occurs  close  to 
compressor  stall  and  unstable  flow  conditions.  In  actual  operation  a  stall 
margin  must  be  provided  to  prevent  the  compressor  from  penetrating  the  stall 
boundary  and  developing  destructive  unsteady  flow  phenomena  such  as  rotating 
stall  and  surge.  This  is  usually  done  by  prescheduling  the  engine  controls. 

When  an  aircraft  has  a  varied  flight  envelope,  the  prescheduling  approach  can 
lead  to  the  requirement  for  a  large  stall  margin  to  keep  the  engine  out  of 
stall  under  all  possible  transient  and  steady  flight  conditions.  This  stall 
margin  represents  a  significant  performance  penalty.  Also,  in  many  instances 
of  engine  failure,  rotating  stall  has  been  identified  as  a  precursor  to  de¬ 
structive  unsteady  flows  in  an  engine.  Furthermore,  blade  fatigue  considera¬ 
tions  will  not  allow  a  compressor  to  operate  for  prolonged  periods  in  a  large 
amplitude  rotating  stall  mode.  Clearly  then,  it  is  desirable  to  develop  methods 
of  estimating  the  stall  boundaries  of  a  compressor  and  if  possible  to  develop 
an  engine  control  system  that  can  sense  incipient  destructive  unsteady  flows 
in  a  compressor  and  take  corrective  action  to  prevent  compressor  stall. 
Recognition  of  these  goals  has  been  the  motivation  for  a  continuing  program 
of  research  that  the  AFAPL  has  sponsored  at  Calspan  since  1962.  This  report 
covers  a  portion  of  the  work  completed  under  the  latest  three  year  research 
program  at  Calspan. 

The  work  accomplished  under  this  portion  of  the  program  was  both 
theoretical  and  experimental  in  nature  and  has  been  aimed  at  obtaining  a  suf¬ 
ficient  understanding  of  the  rotating  stall  phenomena  such  that  its  onset  and 
properties  can  be  predicted  under  a  wide  set  of  circumstances.  The  theoretical 
research  is  reported  in  Section  II  of  this  report  and  covers  the  development 
of  two-dimensional  stability  theories  for  compressible  flows  and  the  progress 
towards  the  development  of  a  three-dimensional  theory  for  incompressible  flow. 
The  experimental  research  is  covered  in  Section  III  of  this  report  and  describes 
a  series  of  three  tests  that  were  performed  in  Calspan/Air  Force  Annular 


1 


Cascade  Facility.  The  first  series  of  tests  was  concerned  with  the  effects 
of  rotor-stator  interference  on  closely  coupled  high  hub-to-tip  ratio  stages. 
The  effects  of  rotor-stator  interference  on  both  the  steady-state  performance 
of  the  blade  rows  and  the  rotating  stall  boundary  for  the  stage  were  determined 
The  second  series  of  tests  was  concerned  with  the  influence  of  a  stationary 
distortion  pattern  on  the  stall  boundary  for  these  stages.  The  third  series 
of  tests  was  concerned  with  obtaining  the  stall  boundary  on  a  low  hub-to-tip 
ratio  rotor.  Section  IV  presents  comparisons  between  the  experimental  results 
and  the  predictions  of  a  two-dimensional  rotating  stall  stability  theory  which 
was  developed  previously.  The  comparisons  include  an  isolated  rotor,  a  rotor- 
stator  stage,  and  the  stage  with  distortion.  The  theoretical  effect  of  chang¬ 
ing  the  axial  separation  between  the  blade  rows  in  a  stage  is  also  discussed. 

A  summary  of  the  program  is  presented  in  Section  V. 


SECTION  II 

THEORETICAL  ROTATING  STALL  RESEARCH 


A.  INTRODUCTION 

The  primary  effects  that  have  been  studied  theoretically  are  the 
stability  of  two-dimensional  compressible  flow  through  a  blade  row  and  the 
stability  of  three-dimensional  incompressible  flow  through  a  blade  row.  The 
formulation  of  a  stability  theory  for  wholly  subsonic  two-dimensional  com¬ 
pressible  flows  given  in  Reference  1  has  been  modified  through  a  revision  of 
the  energy  matching  condition.  Also  some  analysis  of  the  theory  has  been  made 
with  respect  to  the  mechanisms  of  rotating  stall  in  supersonic  flows.  This 
work  is  described  in  Part  B  of  this  section.  Progress  towards  developing  three- 
dimensional  stability  theory  for  incompressible  flow  is  given  in  Part  C. 


B.  TWO-DIMENSIONAL  COMPRESSIBLE  FLOW  STABILITY  THEORY 


A  two-dimensional  stability  theory  for  the  incompressible  flow  through 
a  blade  row  (and  stage)  was  given  in  Ref.  2.  The  theory  employed  a  finite¬ 
thickness  actuator  model  of  the  blade  row  and  required  the  steady  state  loss 
and  turning  performance  of  the  blade  row  as  inputs.  Correlations  of  the  theory 
with  low  speed  data  on  rotating  stall  boundaries  taken  in  the  CALSPAN/USAF 
Annular  Cascade  Facility  have  been  good.  The  theory  has  been  correlated  with 
data  on  isolated  blade  rows  for  sixteen  cases  involving  both  rotors  and  stators, 
and  for  two  different  rotor- stator  stages.  The  theory  generally  predicted 
the  relative  swirl  to  the  blade  row  at  which  rotating  stall  occurred  to  within 
two  percent.  The  basic  theory  was  extended  to  wholly  subsonic  compressible 
flow  in  Ref.  1.  The  compressible  flow  stability  theory  was  shown  to  reduce 
to  the  incompressible  limit  as  but  was  not  correlated  with  data  be¬ 

cause  of  the  lack  of  suitably  detailed  measurements  to  use  as  input  to  the 
theory.  In  the  present  work  the  formulation  of  the  compressible  flow  theory 
given  in  Ref.  1  has  been  modified  in  that  a  revised  energy  matching  condition 
was  developed.  In  addition  some  analysis  of  the  theory  has  been  made  with 
respect  to  the  mechanisms  of  rotating  stall  in  supersonic  flows. 


--  -  ■■  ■■  — -  ,.,M»  ...tV. i  -i  - - 
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A  portion  of  the  analysis  given  in  Ref.  1  will  be  recounted  briefly 
for  the  sake  of  completeness  and  the  ensuing  analysis.  The  blade  row  is  modeled 
by  a  finite  thickness  actuator  shown  in  Figure  1.  The  mean  flows  upstream  and 
downstream  of  the  actuator  have  uniform  swirl.  The  upstream  flow  region  is 
denoted  by  1  and  the  downstream  region  by  2.  Briefly  the  assumptions  inherent 
in  this  type  of  modeling  are: 

1.  The  blade  row  has  high  hub-to-tip  ratio 

2.  The  annular  area  is  constant 

3.  The  flow  disturbances  have  wave  lengths  in  the  circumferential 
direction  (^.-direction)  which  are  large  compared  to  either  the 
blade  spacing  or  chord. 

4.  The  local  performance  of  the  blade  row  is  only  a  function  of  the 
local  inlet  swirl  and  dynamic  pressure  and  may  be  obtained  in  a 
quasi-steady  fashion  from  the  steady  state  performance  data. 


Quantities  in  a  blade  fixed  coordinate  system  are  denoted  with  a 
subscript  o.  The  transformation  between  blade  fixed  and  duct  fixed  coordinates 
is  given  by 


=  * 


where  is  the  blade  velocity  in  the  duct  fixed  system.  This  velocity  is 
given  in  non-dimensional  form  by  jQ.^  *  where  i  denotes  the  flow  region 
and  Ui  is  the  mean  axial  flow  in  region  L  .  The  mean  absolute  swirl  angle  in 
each  region  is  denoted  as  j9^  and  the  tangent  of ^>L  is  defined  as  Si  the  ab¬ 
solute  swirl.  The  relative  swirl  with  respect  to  the  blade  row  is  given  by 
» 

Following  Refs.  1  and  2,  all  of  the  total  flow  quantities,  which  are 
denoted  by  a  hat,  are  decomposed  into  steady  and  unsteady  parts  as 
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Ui=  UL  *  uL(t) 

V/L  s  Wi  +■  ljJ^(t) 

n  -  Pi  -  -fid) 

fi  -  Rl  +  p  Sd 

where  P  and are  the  pressure  and  density,  respectively,  and  the  L  denotes 
the  flow  region.  The  steady  parts  of  each  quantity  are  constants.  The  assump¬ 
tions  on  the  unsteady  flow  are  basically  the  same  as  in  the  incompressible 
case,  i.e.,  the  flow  outside  blade  rows  is  considered  inviscid,  so  that  fric¬ 
tional  and  heat  conduction  effects  are  neglected.  Generally,  however,  the 
flow  will  be  rotational  to  reflect  the  losses  incurred  in  the  blade  row.  The 
unsteady  disturbance  quantities  are  assumed  to  be  much  smaller  than  the  steady 
quantities  so  that  the  equations  of  motion  can  be  linearized  about  the  steady 
flow. 

The  resulting  linearized  equations  of  motion  are  those  of  continuity, 
x-  and  y-  momentum  and  energy;  namely. 


(1) 

(2) 


3x  L  Ri 


=  O 


H.  c  Zx  dp 


=  o 


(3) 

(4) 


where  the  speed  of  sound  is  given  by  «  6  Pi/RL  and  is  constant  in  each 
flow  region;  X  is  the  ratio  of  specific  heats.  The  major  difference  between 
this  problem  and  the  corresponding  one  for  incompressible  flow  is  the  appear¬ 
ance  of  one  more  unknown,  namely  the  density,  and  one  more  governing  equation, 
the  energy  equation.  Eqs.  (1)  to  (4)  can  be  solved  by  assuming  solutions  of 
the  form 
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(5) 

*  «£iC*)e 

(6) 

(7) 

(8) 

Here  C  3  CR+j  C-L  and  n.  is  the  number  of  stall  cells  and  r  is  the  mean  radius 

of  the  blade  row.  The  resulting  disturbance  is  spatially  periodic  in  the  in¬ 
direction  with  period  Ztcn /r  .  The  flow  is  neutrally  stable  if  C ^  -  0, 
unstable  if  O  and  stable  if  ?i>  O  . 


Since  the  system  of  differential  equations  has  constant  coefficients 
the  solutions  for  X-  and  71  i  are  expected  to  be  of  the  form  e 

Using  this  form  in  Eqns.  (5)  through  (8)  and  substituting  into  Eqs.  (1) 
through  (4)  results  in 


Xl  -  DL«3ix  -  At 


oC,i 


RiUL&L-«9l) 


p..  _ 

L  tftt/i  £«t— ‘*0 


*/i.* 


(9) 


7l, 


X 

=  AL  e. 

-b  5Le 

‘-krTnt 

-  E,  S"x 

L 


cii) 

(12) 


where 


°Ln 


+^{/  -  M‘0 + 
I  -  MS 


(13) 


&ZL  = 


THim,  (I  +tJt  >",*)  J 

/  -  MS 


-  --it r(c  "w^) 


(14) 


(IS) 


and 


/m, 


Hi  « 


ix  (c  ^w;) 


*  A  ot- 


3L 


Ux  (axial  Mach  number) 

°-L 

The  A,  ,  3,  ,  D,  and  E,  are  constants  to  be  determined  by  the  match¬ 
ing  conditions.  Normally,  since  there  are  four  simultaneous  differential  equa¬ 
tions  that  we  are  solving  there  would  be  four  values  of  theQ^  obtained. 
However  for  this  particular  system  0(3^is  a  double  root.  Nevertheless  four 
linearly  independent  solutions  can  be  obtained  for  this  case  using  the  three 
■ 
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The  terms  of  the  disturbance  quantities  with  exponential  factors 
06c  and OCZi_  correspond  to  irrotational  disturbances.  Those  terms  of  U^. and 
with  exponential  factors  correspond  to  vorticity  waves  ccnvected  with 
the  mean  flow  speed.  The  El  term  of  the  density  disturbance  corresponds  to  a 
convected  entropy  wave  as  will  be  seen  by  the  following.  For  a  perfect  gas 
the  entropy, ,  in  flow  region  i  is  given  by  the  following  expression 


_R 

M 


JU 


(16) 


where  -4^  is  a  reference  value  and  R  is  the  gas  constant  for  the  particular  gas 
under  consideration.  Using  the  linearizing  assumptions  this  becomes 


aa  -h 


JB-{ 

r-i  c 


Mu  Pi  + 

W 


—  (fL  -  j 


(17) 


The  onlv  non-constant  terms  in 


-  CL 


in  this  expression  are  theytj^and 
oC,.  X-f 

if>L  »  -£) 


terms  and 


(18) 


At  this  point  it  should  be  noted  that  disturbance  solutions  corre¬ 
sponding  to  Eqs.  (9)  through  (12)  have  been  given  by  Smith.3  The  solutions 
presented  in  Ref.  3  are,  however,  in  error  by  the  omission  of  the  term  in 
the  density  perturbation  (the  entropy  wave  term).  However  the  entropy  and  the 
vorticity  are  related  by  Crocco's  theorem  which  is  given  in  Ref.  4  as. 

+  V  xg  =  g^Ji*  (19) 

where  7"is  the  temperature,  V  is  the  velocity  vector,^  is  the  vorticity 
vector  and-^,d  is  the  total  enthalpy.  Hence,  it  is  not  consistent  to  retain 
the  vorticity  disturbance  terms  without  the  entropy  disturbance  terms  in  the 
disturbance  solutions. 

Thus  far  the  solution  for  the  disturbance  velocities  contain  eight 
undetermined  constants  namely  the  Ai  ,  B ^  ,  0^  and  E{_ .  For  wholly  subsonic 
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flow  these  constants  are  determined  by  requiring  bounded  flow  far  upstream 
and  downstream,  requiring  irrotational  flow  upstream  of  the  blade  row,  and  by 
applying  matching  conditions  across  the  blade  row.  The  requirement  for  bounded 
flow  far  upstream  and  downstream  results  in  B^=  O  .  The  requirement  for 
irrotational  flow  upstream  of  the  blade  row  results  inD(=£,=0.  The  remaining 
four  constants  are  determined  by  applying  matching  conditions  across  the  actu¬ 
ator  which  are  derived  from  the  continuity  equation,  the  axial  momentum  equa¬ 
tion  (with  losses),  the  energy  equation  and  the  blade  row  steady-state  turning 
performance. 


The  matching  conditions  are  derived  using  the  notation  and  symbols 
shown  in  Figure  1  and  are  generally  written  in  the  blade  fixed  coordinate 
system.  The  conservation  of  mass  flow  across  the  actuator  gives  the  mean  flow 
quantities  Ut  R,=  (-4R2and  for  the  unsteady  quantities 

* 


J 


U,^Ol  -t  U.,R ,  -  Jyoa  -  uzRz  +  cos  S'  S  =  O 


(20) 


The  quantities  with  a  subscript  1  are  evaluated  at  X0(  =  O  and  l^QI  ;  the  quan¬ 
tities  with  a  subscript  2  are  evaluated  at  X0;t  =  Jl  cos  S  and  =  <^o,+-(?scn  S 
The  integral  term  in  Eq.  (20)  is  evaluated  by  assuming  a  linear  variation  in 
the  integrand  between  O  and  . 


It  is  assumed  that  the  turning  performance  of  the  blade  row  may  be 
expressed  as  . 

J,  -  G(X)  (2D 

where  G  represents  the  steady  state  turning  of  the  blade  row  functionally. 
Linearizing  this  relation  about  its  steady-state  value  gives 


(22) 


wp.ere 
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The  vector  form  of  the  momentum  equation  may  be  put  in  the  following 


form 

TT 


v^o)  ■  dt  +  ^  f-cL7 


(23) 


where  clr  is  the  differential  of  the  position  vector  to  points  on  the  contour 
C  in  Figure  1.  The  line  integrals  are  taken  around  this  contour.  The  fric¬ 
tion  force  vector,  K  ,  is  most  easily  defined  in  Cortesian  Tensor  notation. 

Let  X,=  X  and  X2-  tc  then  the  ith  component  of  F  is  3  ^  where  the  summation 

0  ~  r  °xj 

convention  should  be  used  and  1  Cj  is  the  viscous  stress  acting  in  Xj  direction 
on  planes  perpendicular  to  .  For  laminar  flows  the^ij  may  be  computed 
using  the  Newtonian  assumption  and  Stokes  hypothesis.  For  turbulent  flows 
the  viscous  stress  are  usually  given  empirically.  Eq.  (23)  may  be  evaluated 
as  in  Ref.  2  to  give  ^ 

(y°rw°z)~s&c S  xCo  5 "fir  z 


-  U, '  [x+  4(1  +  J?)  x]  +  [xij-  xO  X' ] 


(24) 


T 


=  o 


where  X  is  a  loss  coefficient  for  steady  flow  defined  by. 

n  .  — CPT.,  -  FtO 


where  Pra.  is  the  steady-state  total  pressure  relative  to  the  blade  row  in 
flow  region  i  .  Functionally,  it  is  assumed  that 


so  that 


x  ■  xQS,) 


x 


ZJ, 


(25) 


The  remaining  matching  condition  comes  from  the  energy  equation. 

An  energy  matching  condition  was  derived  in  Ref.  1,  however,  this  matching  con¬ 
dition  required  the  evaluation  of  the  integral  of  a  dissipation  function  through 
the  blade  row.  This  was  accomplished  by  assuming  linear  variations  of  the 
quantities  to  be  integrated  based  upon  their  values  immediately  upstream  and 
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downstream  of  the  actuator.  This  assumption  has  also  been  used  to  evaluate 
integrals  appearing  in  the  continuity  and  momentum  matching  conditions,  but 
is  not  appropriate  for  evaluating  the  integral  of  the  dissipation  function. 

The  dissipation  function  involves  the  product  of  the  viscosity  and  squares 
of  derivatives  of  the  velocity  components.  In  the  actual  flow  through  the 
blade  row  these  derivatives  may  be  large  thus  estimating  these  derivatives 
from  the  model  of  the  upstream  and  downstream  flows  which  are  composed  of 
constant  values  for  the  mean  flow  velocities  plus  small  unsteady  components 
is  inaccurate.  Also  the  form  of  the  dissipation  function  used  is  only  strictly 
valid  for  laminar  flows.  A  form  of  the  energy  equation  has  been  found  which 
avoids  evaluation  of  the  dissipation  function  and  will  now  be  derived. 

The  starting  point  is  the  basic  energy  conservation  equation  given 
in  Ref.  5.  Neglecting  heat  transfer  by  conduction  this  equation  is 


Wy. 


(26) 


and  is  applied  to  an  arbitrary  control  volume,  V0  ,  bound  by  the  surface 
The  Cartesian  tensor  notation  and  summation  convention  are  used  here  again. 

The  integrals  appearing  on  the  left  hand  side  of  the  equations  are  surface 
integrals  where  is  the  ith  component  of  the  velocity  vector's/  ,  is  the 
ith  component  of  the  outward  normal Y\  at  the  bounding  surface.  The  stresses 
at  the  control  volume  surface  are  i>ij  and  are  composed  of  the  pressure  forces  and 
viscous  stresses  such  that 


when- 


LJ 

is  the  static  pressure;  Ti 


Tj  7 o  $Lj 

are  the  viscous  stresses  previously  defined 


-  —  - - - - -  ,tJ 

and  S;;  is  the  Kronecker  delta  function  having  the  property 

J  Sij-l  -i-J 


=  O  S.  \  j 


£".£  is  the  sum  of  the  internal  energy,  e.  ,  and  kinetic  energy  per  unit  mass 
of  fluid  so  that 
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and  e*CvT for  a  thermally  perfect  gas  where  T~ is  the  temperature  and  Cv  is 

^  d. 

the  specific  heat  at  constant  volume,  "^fis  the  rate  of  heat  produced  per  unit 
volume  in  M>  and  is  due  to  friction  in  this  case. 

Basically  Eq.  (26)  states  that  the  rate  of  work  performed  on  the 
control  volume  by  the  surface  stresses  minus  the  rate  at  which  energy  is 
transported  out  of  the  volume  is  equal  to  the  time  rate  of  change  of  energy 
inside  the  volume  minus  the  heat  produced  in  the  volume.  To  apply  Eq.  (26) 
to  the  present  problem,  choose  a  control  volume  (area  for  the  two-dimensional 
case)  bounded  by  the  curve@in  Figure  1.  In  the  figure,  ©is  drawn  so  that 
B-C  and  A-D  correspond  to  blade  surfaces.  This  is  not  necessary  but  these 
two  sides  of  the  control  volume  should  be  chosen  such  that  they  are  parallel 
to  the  blade  chord  lines.  The  assumption  will  be  made  that  flow  in  the  passage 
is  parallel  to  the  blade  chords  so  that  convection  of  energy  out  these  sides 
of  the  control  volume  can  be  neglected. 


We  now  use  an  argument  given  by  Vavra  that  the  frictional  heat  is 
immediately  and  equally  transmitted  to  all  points  in  a  fluid  particle  as  it 
moves  along  its  path.  Further  it  is  assumed  that  the  frictional  heating  is 
equivalent  to  the  work  of  the  frictional  forces  acting  on  the  surface  of  the 
particle.  Under  these  assumptions  for  the  control  volume  chosen  the  viscous 
portion  of  the  surface  stress  work  integral  (first  integral  on  left  hand  side 
of  Eq.  26)  is  cancelled  by  the  heat  production  integral  (last  integral  on  right 
hand^side  of  Eq.^26).  Then  Eq.  (26)  may  be^simplified  to 


•J 


(27) 


i  U'£d}  -J fz  U* 

=  C°SS 

where  are  a  sheared  coordinate  system  parallel  to  the  blade  chord  and 

cascade  axis  respectively.  Since  £  is  arbitrary 

^  V*'  &  £t,  7 =  cos  S  J  £t)  ds  (28) 


/ 


: 
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After  applying  the  linearizing  assumptions  and  use  of  the  mean  flow  continuity 
equation  this  equation  can  be  split  into  a  steady  part 


yrr  ?i ;  ^[u,s(w,-  wt)* 


(29) 


which  essentially  states  that  the  mean  total  enthalpy  is  conserved  and  an 
unsteady  part 


U>t' Tfr}+«.{R.uNE.R,]+^  R(u((w,-wb) 

'  hi  h  -v]  ”  ^  “  W2  -wb) 

•"  ^ coi  s  +/i'(5<,'c^S.)  +^'R,U>  +  ^-Wb) 


+  A. 

r-i 


V*  (£-&■) + U’R*U* + 


f?2  (w,  -  w  b  )T 


(30) 


Here  a  dot  over  a  quantity  signifies  its  time  derivative  and  the  integral 
appearing  on  the  right  hand  side  of  Eq.  (28)  has  been  evaluated  by  assuming 
a  linear  variation  of  the  integrand.  Also  we  have  let 

having  the  same  value  for  both  i  -  1  and  2  because  of  Eq.  (29).  (This  is  the 
mean  total  enthalpy  of  the  flow.) 


It  will  be  noted,  by  inspection  of  Eq.  (29)  and  (30),  that  the  vis¬ 
cosity  does  not  explicitly  appear  in  the  energy  matching  conditions.  This 
appears  to  be  consistent  with  the  energy  equation  used  for  one-dimensional 
flow  in  pipes  with  friction,  (see  Ref.  7). 

Eqs.  (20),  (22),  (24)  and  (30)  supply  the  required  relation  to  de¬ 
termine  the  unknown  constants.  Substituting  Eqs.  (9)  through  (12)  in  these 
equations  results  in  a  system  of  homogeneous  equations  of  the  form 
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[A]  [B,  1  =  O 

Da  ( 

Ai  1 

e :2 

where  the  A  matrix  has  elements  A-j  .  The  characteristic  equation  of  the 
svstem  which  gives  the  allowable  values  of  C  is  then 

det  [A]  =  O  1 

The  elements  of  the  A  matrix  are: 

A,  -  *  ££*< 

R  U* 

A  -  '  ' 

***  uz 

A  »  n*  —  -  i  a  (— )  7-^  cos  ^  -  ~u~f~ 
z  uz  f  z  [*z'  r  <44 

A*  ■  Mz  (,  +  *T  7£  CO50 

^2/  *  0  -  M  j  "  X?  J*} 

(£-4) 

AZ3  *  7^  {_4J3  0-^j} 


-  0  -  hi)  £•*■  +  xj,  +  ~0 1 
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5  (U  *  Tf)  *  j  4  ft-  ™ a)  u, 

-  j  s  f-%  ■»  — l 

2  i  »*  ] 


A«  ■  -^C4  -c.% 


A-,3  3  -  U> 


^3  (Uz'-Vt) 

1  +J 

,4i 

(/-tf 

>4 

I 

*  ♦.»' 

—  1 
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5- 
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1 
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where 


I,  •  4 

I,  •  jf,  *  [l 

I,  •  j.  M,' M*  0  *  O] 
i4  •  ??,-  [i-«,*C'  +  £)|a 

A  *  -T7T-  *  ■&. 


M,  -  U'A 

H*  *  V*. 

As  may  be  seen  from  inspecting  the  Aij  elements  of  the  A  matrix 
Eq.  (32)  is  not  a  simple  polynomial,  as  in  the  incompressible  flow  case, 
and  must  be  solved  numerically.  This  has  not  been  done  yet.  However,  Eq. 

(32)  reduces  to  the  incompressible  limit  of  Ref.  1  when CL^-*-co 
and  the  numerical  results  should  be  the  same  as  those  presented  in  Ref.  1  for 
incompressible  flow. 


As  mentioned  previously  the  foregoing  analysis  is  limited  to 
wholly  subsonic  flows.  This  limitation  results  because  of  the  steady-state 
blade  row  performance  assumed.  The  basic  solutions,  Eqs.  (9)  through  (12) 
for  the  disturbance  quantities  are  valid  for  supersonic  flows.  Examining  the 
exponential  factors  c^ti  and  OLzi  (Eqs.  13  and  14)  appearing  in  the  expressions 
for  the  velocity  disturbances  for  the  case  of  neutrally  stable  flows  ( C x  =  0^ 
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it  can  be  seen  that  the  variation  with  X  changes  character  when 

U2  +(v/, -V^,)2  >  a.  2  (33) 

Here  is  the  propagation  velocity  of  the  wave  in  the  azimuthal  direction. 

When  this  inequality  occurs  the  radical  terms  of  Eqs.  13  and  14 
change  from  real  to  imaginary  and  the  variation  with  X  changes  from  damped 
oscillatory  (complex  valued cXlJ  )  to  undamped  oscillatory  (purely  imaginary 
valued  oCcj  ).  At  this  point  the  disturbances  quantities  correspond  to  sound 
waves  that  are  propagated  to  upstream  infinity.  The  selection  of  which  values 
of  thecX,;:  apply  in  each  flow  region  is  then  determined  by  the  requirement  that 

V 

the  waves  emanate  at  the  blade  row  and  travel  away  from  it  in  both  the  upstream 
and  downstream  direction.  Put  another  way,  no  waves  should  emanate  from  either 
upstream  or  downstream  infinity.  Assuming  that  W-Vp>0,  the  OiZ[  root  is  used 
in  the  upstream  flow  and  the  <XI2 root  is  used  in  the  downstream  flow  as  before. 


The  condition  expressed  by  relation  (33)  occurs  when  the  inlet  flow 
becomes  supersonic  relative  to  the  disturbance  wave.  Since  Wb  is  negative  for 
a  compressor,  in  the  sign  convention  being  used,  and  since  is  usually  about 
one  half  ,  relation  (33)  will  generally  occur  at  some  point  after  the  inlet 
flow  has  become  supersonic  relative  to  the  blades. 

There  are  several  flow  regimes  to  consider  for  nominally  supersonic 
flows  dependent  upon  whether  the  axial  flow  is  supersonic  or  not.  Further  the 
cases  with  subsonic  axial  flows  have  different  characteristics  dependent  upon 
whether  both  the  inlet  and  outlet  flows  are  supersonic  relative  to  the  rotor. 
The  steady-state  characteristics  of  cascades  operating  in  these  different  flow 
regimes  is  discussed  in  Ref.  8.  In  the  present  discussion  only  the  case  with 
supersonic  axial  flow  and  the  case  with  subsonic  axial  flow  where  both  the 
inlet  and  outlet  conditions  are  supersonic  relative  to  the  rotor  will  be  con¬ 
sidered. 


First  the  case  with  subsonic  axial  flow  will  be  discussed.  The 
steady  operational  characteristics  of  the  cascade  have  two  properties  which  are 
very  important  to  this  case.  First  of  all,  a  supersonic  cascade  with  subsonic 


17 


axial  flow  can  only  operate  in  the  steady  state  at  a  single  inlet  angle. 

This  condition  is  often  referred  to  as  "unique  incidence".  Secondly  the  back 
pressure  (flow  angle)  far  downstream  of  the  cascade  is  not  a  unique  function 
of  mass  flow  and  inlet  conditions  to  the  cascade.  A  range  of  back  pressures 
is  possible  without  unstarting  the  cascade;  therefore,  the  back  pressure  must 
be  prescribed  as  part  of  the  boundary  conditions  for  the  problem.  In  addition 
there  is  only  one-way  communications  for  the  unsteady  disturbances  between 
upstream  and  downstream  flow  fields  (assuming  that  the  disturbances  are  small 
enough  not  to  unstart  a  blade  passage) .  Nothing  that  happens  in  the  downstream 
flow  can  influence  the  upstream  flow.  Propagation  of  a  disturbance  which 
starts  at  the  trailing  edge,  upstream  through  the  blade  passage,  is  blocked  by 
the  blade  surfaces  and  reflected  back  into  the  downstream  flow.  Therefore 
the  upstream  disturbance  field  may  be  neglected  as  far  as  stability  considera¬ 
tions  are  concerned  and  there  are  only  three  constants  (A^,  and  E_J  to  be 
determined.  This  requires  the  following  revisions  to  the  required  matching 
conditions.  First  the  turning  relation,  Eq.  (22),  is  dropped.  Secondly  the 
loss  relation,  Eq.  (25),  should  be  revised  to  indicate  that  now  the  losses 
are  a  function  of  the  downstream  flow  angle  (or  static  pressure)  and  this  re¬ 
vision  incorporated  in  the  vorticity  compatability  matching  condition,  Eq.  (24) 
Thirdly  the  three  matching  conditions  used  should  be  revised  to  reflect  the 
fact  that  all  the  upstream  disturbances  are  now  zero.  The  characteristic 
equation  for  this  case  can  be  obtained  quite  simply  from  Eq.  31  by  deleting 
the  first  column  and  second  row  of  the  A  matrix.  In  addition  the  third  row 
must  be  revised  because  of  the  new  loss  relationship.  The  resulting  equation 
is  however,  again,  not  a  simple  polynomial  and  must  be  solved  numerically. 

This  has  not  been  done  yet. 

For  the  case  of  supersonic  axial  flow  the  unique  incidence  does  not 
apply  and  the  outlet  conditions  are  uniquely  determined  by  the  inlet  conditions 
However,  there  can  be  no  disturbances  ahead  of  the  cascade  so  that,  again, 
there  are  only  three  undetermined  constants  to  be  determined  in  the  disturbance 
solutions;  namely  A.,,  D 2  and  E2.  At  this  stage  is  not  clear  which  matching 
condition  should  be  dropped.  It  appears  that  this  question  will  have  to  be 
settled  by  correlation  with  experimental  data. 
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These  considerations  of  the  different  types  of  compressible  flows 
show  that  the  two-dimensional  analysis  should  be  limited  to  cases  where  the 
same  type  of  flow  prevails  along  the  entire  span  of  the  blade.  Moreover  it 
appears  that  the  most  significant  blade  row  performance  parameters  that  in¬ 
fluences  the  inception  of  rotating  stall  must  be  substantially  different  in 
the  supersonic  cases.  It  was  concluded  in  Ref.  2  that  the  slope  of  blade  row 
losses  with  inlet  flow  angle  was  the  most  important  performance  parameter 
with  respect  to  the  stall  boundaries  for  incompressible  flows  (this  also 
appears  true  for  wholly  subsonic  compressible  flows).  For  the  case  of  super¬ 
sonic  flow  with  subsonic  axial  flow  the  blade  row  losses  will  be  dependent 
upon  the  back  pressure  rather  than  the  inlet  flow  angle  which  suggests  that 
the  slope  of  the  blade  row  losses  with  back  pressure  may  be  the  controlling 
factor  for  this  case.  The  case  with  supersonic  axial  flow  remains  a  subject 
for  further  analysis. 

C.  THREE-DIMENSIONAL  INCOMPRESSIBLE  FLOW  STABILITY  THEORY 

There  are  basically  two  ways  to  study  the  stability  of  a  system. 

The  first  way  is  to  determine  under  what  conditions  small  disturbances  with 
an  oscillatory  time  behavior  are  either  damped  or  amplified  with  time.  This 
is  referred  to  as  a  stability  theory  approach  and  was  used  in  the  two-dimen¬ 
sional  case.  The  second  method  is  to  disturb  the  system  and  examine  the  time 
history  of  the  response  of  the  system  to  determine  if  the  response  becomes 
unbounded  with  time.  The  second  method  is  somewhat  more  general  in  that  it  is 
not  necessarily  limited  to  linear  systems.  However,  it  involves  actually 
calculating  the  response  of  the  system  for  a  variety  of  input  disturbances. 

If  this  is  done  numerically,  care  must  be  taken  to  ensure  that  divergence  of 
the  response  is  actually  due  to  the  physics  of  the  problem  and  not  due  to 
instabilities  of  the  numerical  algorithm  used.  It  is  therefore,  usually, 
better  to  use  a  stability  analysis  if  it  is  tractable. 

In  the  two-dimensional  case  a  stability  theory  has  proved  quite 
useful  in  understanding  the  mechanisms  of  rotating  stall  and  predicting 
inception  conditions  for  blade  rows  with  high  hub-to-tip  ratios.  However, 
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most  modern  compressors  have  initial  stages  with  relatively  low  hub-to-tip 
ratios  and  it  is  thought  that  three-dimensional  effects  may  be  important  for 
these  cases.  The  two-dimensional  theory  for  an  isolated  blade  row  provided 
concise  results  in  that  for  the  incompressible  case  a  closed  form  expression 
was  obtained  for  the  damping  factor  and  propagation  velocity  of  admissible 
disturbance  solutions  (Eqs.  30  Ref.  2).  In  addition  the  conditions  that  apply 
at  the  neutral  stability  boundary  could  be  concisely  stated  (Eq.  32  Ref.  2). 
Therefore  the  stability  boundary  for  a  given  blade  row  could  be  calculated 
with  relatively  little  numerical  effort.  The  case  for  two  blade  rows  was  only 
slightly  more  complicated  so  that  it  would  be  practical  to  apply  this  theory 
to  several  stages  if  the  appropriate  blade  row  loss  and  turning  performance 
were  available.  The  conciseness  of  the  results  for  the  two-dimensional  case 
also  promoted  a  physical  interpretation  and  understanding  of  the  results. 

This  situation  prompted  the  formulation  of  a  three-dimensional  stability  anal¬ 
ysis  using  the  same  overall  approach  as  in  the  two-dimensional  case.  The 
simplified  model  for  the  mean  flow  developed  in  Ref.  1  was  used.  However 
applying  the  overall  methodology  of  the  two-dimensional  problem  to  the  three- 
dimensional  case  resulted  in  a  much  more  complicated  problem  that  did  not  allow 
a  concise  statement  of  results.  Further,  sufficient  numerical  effort  would 
have  been  required  such  that  physical  insight  to  the  problem  would  be  lost. 
Therefore  the  implementation  of  the  three  dimensional  theory  was  not  completed. 
However  the  theory  has  been  developed  to  the  point  where  it  could  be  imple¬ 
mented  and  this  development  is  presented  in  the  following  paragraphs. 

General  Analysis 

As  in  the  two-dimensional  case  the  basic  plan  of  attack  was  to  deter¬ 
mine  basic  solutions  for  the  perturbation  velocities  and  pressures  both  up¬ 
stream  and  downstream  of  the  blade  row.  It  was  expected  that  these  disturbance 
solutions  would  contain  unknown  constants  that  would  be  determined  by  the 
matching  conditions  across  the  blade  row.  These  matching  conditions  would 
result  in  a  characteristic  equation  for  the  damping  factors  and  propagation 
velocities  for  admissible  solutions.  The  major  difficulty  arose  when  the 
downstream  disturbance  velocities  could  not  be  determined  in  analytical  form. 
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Coupling  a  numerical  solution  for  the  downstream  disturbance  velocities  with 
the  characteristic  equation  obtained  from  the  matching  conditions  leads  to 
the  requirement  for  substantial  numerical  effort  to  implement  the  analysis. 

The  implementation  was  not  completed  but  was  carried  far  enough  that  the  de¬ 
tailed  mechanics  of  obtaining  a  solution  are  relatively  clear.  The  completed 
portion  of  the  analysis  will  be  given  along  with  the  intended  method  of  solu¬ 
tion. 

The  simplified  actuator  sheet  model  developed  in  Ref.  1  was  used 
as  the  model  for  the  basic  steady  flow.  Basically  the  flow  considered  is  that 
in  a  constant  area  annulus.  The  blade  row  is  modeled  by  an  actuator  sheet. 

The  inner  radius  of  the  annulus  is  denoted  by  rH  and  the  outer  radius  by  vT  . 
The  following  sketch  shows  the  coordinate  system  and  geometry. 


The  mean  flow  is  axially  symmetric.  The  axial,  radial  and  azimuthal 
velocity  components  are  designated  by  U,  V  and W  respectively.  From  Ref.  1 
we  have  that  a  stream  function  tA/  exists  such  that 

v---A^8  i  U.  -1-  w 

and  the  approximate  solutions  for  and  the  other  flow  quantities  are  given 
by 


V'-U.-J-‘  +  a,(»)('K..Kirtr.) 

W-  —  ip  +  -2s. 

v  ~  r 

H=  A  ^  +  D, 


(34) 

(35) 

(36) 
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where  H  is  the  total  pressure  divided  by  the  density,  U®  is  the  undisturbed 
axial  velocity  far  upstream  of  the  blade  row,  and 


Ko  -  Vtp 


The  function 

a2  (2)  has  two  forms  dependent  upon  Z  ?•  O  • 

Upstream 

when 

z<  0 

H.* 

e 

(37) 

and  for 

z  >  o 

where 

* 

—  (M  -r*n  -  Mi  <2  ) 

(38) 

( 

'6K.X.  £ 

M  = 

(w*-  8*)* 

However,  in  the  ensuing  development  we  will  not  distinguish  between  the  two 
forms  of  cl2  until  necessary.  Also  A,  B,  D1  and  are  constants.  Upstream 
of  the  blade  row  A,  B  and  D,  are  zero.  Downstream  of  the  blade  row  they  may 
be  determined  by  using  two-dimensional  cascade  data  in  a  strip  theory  fashion, 
as  in  Ref.  1,  or  they  may  be  determined  to  give  a  best  fit  to  experimental  ve¬ 
locity  distributions  if  they  are  known.  In  any  event,  they  are  considered 
as  known  quantities  for  the  present  analysis. 

The  total  velocity  components  and  pressure,  which  are  denoted  by 
a  hat  are  decomposed  into  steady  and  unsteady  parts  as  follows 


UfoO  +  U.  (a,r)  e 


(39a) 
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(59b) 


w  \/  r  \  ~r  x 
V  -  V  (*,r)  +  v  (s.'O  ft 

w  *  W(k,r)  +  w  (e,v)  e  (39c) 

fc,Y)  +  (*,r}  4  "*  n  ^  (39d) 


The  major  difference  between  these  forms  and  those  used  in  the  two- 
dimensional  case  is  that  the  steady  parts  are  no  longer  constants  but  are  de¬ 
termined  from  Eqs.  (34)  through  (36).  Assuming  that  the  steady  part  of  each 
quantity  is  much  larger  than  the  unsteady  part,  the  momentum  and  continuity 
equations  may  be  linearized  into  the  following: 


+  :.llf  (40) 

'  A  v  A  v  V  '  "*  '  —*  —  — 


dy  Jr  ' J  '  r  '  1  "  J2  ’  " 

jc^+V|2  +  ~4^+jnJ^Lw+u4^  +IX^+v 

ay  ay  b  r  h-a 


/° 


W 

V 


..li£ 


-j'-— 


I  ^  /"  *-v-\  ,  /*v  ^  _ 

~  +j  —  w  Ti  =  o 


(41) 


(42) 

(43) 


Eqs.  (40) .  (41)  and  (42)  are  the  v  ,  0  and  2  momentum  equations  respectively 
and  Eq.  (43)  is  the  continuity  equation. 

Based  upon  the  approximate  techniques  used  in  Ref.  1  solutions  for 
the  disturbance  quantities  are  sought  in  the  form 

ttx  ii^fi+U,C*)  (44.) 
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v  *  +  V,  (*)  +  vz  (2)  r 

w  =  ^-+W,60  -v  w2  (z)  r 

y* 

(*)+?»,  (Z-)r  +-f  2  60Y‘ 


Requiring  no  flow  through  the  annulus  walls  and  the  continuity  equa¬ 
tion,  Eq.  (43)  to  be  satisfied  for  all  r  results  in. 


K,  4-r) 

(45) 

1  +-  [2V,  +  U./  ]| 

(46) 

where  primes  denote  differentiation  with  respect  to  2  .  Thus  there  are  six 
remaining  functions  of  Z  to  be  determined;  namely  U.a  ,  (i.,  >  ,  -Jo<>  >  and 

^a-2  .  The  six  required  equations  to  determine  these  functions  will  be  deter¬ 
mined  in  different  fashions  for  the  upstream  and  downstream  flows. 

Upstream  Flow  Disturbances 

The  upstream  flow  is  much  easier  to  analyze  since  it  is  irrotational 
and  the  velocity  components  may  be  solved  for  independently  of  the  pressure. 
Let  the  ith  component  of  vorticity  be  denoted  by-£2.cand  let 

^  jCct*„e) 

=  ±lL  e 

then 

-£*-  v(£C<-w)-j"v)  («*> 

.  j  _  4*1  (47b) 

^  v  'J  v  £  2. 
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(47c) 


^  v 

c)-Z. 


^  r 


Now  the  equations  required  to  solve  for  (J-ot  U-,  and  Vz  are  obtained 
by  substituting  Eqs;  (44a).  (45)  and  (46)  into  Eqs.  (47)  and  requiring 

J*  il  '  dr  =  0 

This  leads  to  the  following  system  of  ordinary  differential  equations  for 
the  unknown  function 


4  v'  +  is.  u.  r  0 

O 


where 


(48a) 

(48b) 

(48c) 


vx  (2  4la  -  A  O  +  A  +  ZA0U'>  s  0 

u."  +  \  ■  0 


■4®  *  '^v  (rH  ) 

^  *  (  *T  ~  r«  ) 

^,o  •  iW-O 

Since  this  system  has  constant  coefficients  the  solution  will  be  of  the  form 


U.  =  C.c 


»■  C,  C 
V2  = 


A.Z 


Substitution  of  these  forms  in  Eqs.  (48)  gives  the  following  equation  for  X 


*  "  ^  <>*i  +  <*,o)  +  0,o  ~  *Soa»,)  -  0  W 

where 


2 

-  n. 


^6 
■A  to 


The  Teal  positive  roots  of  Eq.  (49)  provide  the  values  of  X  of  interest. 
Once  these  have  been  determined  Eqs.  (48)  provide  two  relations  involving  c  , 

Cj  and  c2>  namely 

and 


Thus  for  each  value  of  X  there  is  one  undetermined  constant  associated 
with  the  upstream  disturbance  velocities  that  remains  to  be  determined  by 
the  matching  conditions.  The  ^  may  be  determined  by  substituting  the  above 
solutions  for  U.„  ,  U- ,  and  V2  into  Eqs.  (40),  (41)  and  (42).  At  this  point 
we  note  that  the  upstream  disturbance  velocities  are  independent  of  the  mean 
flow.  This  is  a  consequence  of  the  irrotationality  of  the  upstream  flow. 

This  fact  also  allows  us  to  construct  an  exact  solution  for  the  upstream-flow 
disturbances  that  may  be  compared  to  the  approximate  solutions  just  found. 

Since  the  upstream  flow  is  irrotational,  a  velocity  potential  (^exists  such 
that 
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where 


/ 


§ m  0C*,r)e 


j (c"t  +  n  ©") 


a=  0Z  ,7:0, 


V- >  WSJ 


n  0 


and  the  continuity  equation  (Eq.  (43))  becomes 

+<A-r  *  -T0V-  -^<0=0 
Then  0  has  a  Fourier  Bessel  expansion  given  by 

-  |0  [v  3.  (<*;  -£)  +  52.  •  y„  («(i)]  e“‘ *  (=0) 

— .  V  "tj( 

Here,  J*  and  /„  are  the  *\  order  Bessel  function  of  the  first  and  second 
kind,  respectively.  The  OCL  are  the  roots  of 

te;K'C«iR)  -  j„' <*)/„'(«;*?)  =  0  ‘51> 

The  primes  denote  differentiation  with  respect  to  r  ,  and  R  is  the  hub-to-tip 
ratio,  Y\i/y~r.  The  constants  S,^  and -S^  satisfy  the  following  relation 


£lL 

s. 


2i 


y »(cclr) 


(52) 


These  constants  are  determined  for  each  cCL  by  prescribing  appropri¬ 
ate  boundary  conditions  at  the  blade  row  (2=0)  which  result  from  matching 
Cl  and  v  with  the  downstream  solution  and  by  prescribing  a  turning  relation 
through  the  blade  row. 


Eq.  51  results  from  the  condition  of  no  flow  through  the  duct  walls 
and  the  Oij_  which  satisfy  this  equation  must  be  determined  numerically.  Tables 
of  these  roots  are  given  by  Bauer  for  0£  V>^/(  and  0  9  . 

A  comparison  of  the  exponential  decay  factors  for  R  =  .8  (correspond¬ 
ing  to  the  Calspan/AF  annular  cascade  tests)  will  now  be  made.  Eq.  49  yields 
two  solutions  for  the  decay  factor  of  the  approximate  solution  which  corre¬ 
sponds  to  the  first  two  terms  in  the  Fourier- Bessel  expansion  of  the  exact 
solution.  For  K  =  1  the  value  of  0do  given  in  Ref.  8  is  1.11337  and  that  for 
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OC,  is  15.77771  indicating  that  the  subsequent  terms  in  the  exact  solution  decay 
very  rapidly  as  distance  from  the  blade  row  is  increased.  The  agreement  between 
the  decay  factors  of  the  exact  and  approximate  solutions  for  the  first  term 
is  excellent  as  shown  in  the  following  table. 


COMPARISON  OF  EXPONENTIAL  DECAY  FACTOR  IN 
UPSTREAM  DISTURBANCE  SOLUTIONS 


1st  Term  in  Exact 
Bessel  Function  Expansion 


Approximate 
Velocity  Components 

~K'rT 


n  =  1  1.11337 
n  =  2  2.22646 
n  =  5  5.56194 
n  =  10  11.07246 


1.11337 

2.22652 

5.56233 

11.09664 


The  axial  decay  factors  for  the  second  term  in  the  approximate  solu¬ 
tion  are  about  ten  percent  high.  However,  overall,  this  type  of  agreement 
between  the  approximate  and  exact  solution  indicate  that  the  approximate  method 
should  be  sufficiently  accurate. 


Downstream  Flow  Disturbances 


For  the  downstream  flow  the  six  unknown  functions  Ua ,  LL,  ,  V*  , 

,  -jo,  and are  determined  by  substituting  Eqs.  (44a),  (45)  and  (46)  into 
Eqs .  (40)  through  (49),  integrating  between  and  ,  and  taking  the  first 
integral  moment  of  these  equations  between  and  r*T  .  This  leads  to  a  system 
of  ordinary  differential  equations  of  the  form 

17/  V2  +  fi!z  Vk  +  /T3  ua  *  /"i4  U-o  t  As  <4-0 
'Tnu,  +rr7  u,"  +  r?8  u"  +  +Jtii  -fa 
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(53) 


+  JL^fSAi<?:+A.LSfz*M(.f'l  -  O 

■for  L  ~  /,  2  j . .  £ 

The  detailed  expressions  for  the  Tjj  an  dhj  are  given  in  Appendix  A.  The 
rjj  are  functions  of  Z  because  &4is  a  function  of  Z  .  The  -/Lij  are  constant 
terms.  This  allows  theyio^  to  be  eliminated  from  the  system  of  equations  ar¬ 
riving  at  the  somewhat  simplified  form 

fi,\ i  +  Hi(  v/  4  P^,  y/  *Futl.*Hi3tu/  *  Pi2  a/' 

+  +  +  Hi3U,\FizU,%Qi:ilL*'  =  O  (54) 

for  i  r/,2,3 

The  FV-  ,  )4ij  ,  PtJ  and  Q  ij  are  functions  of  2  given  in  the 
Appendix  A.  It  will  be  noted  here  that  the  functional  dependence  upon  z  is 


of  a  special,  form  displayed  as  follows: 

F‘J 

(55a) 

(55b) 

Pij  =  pLj  +A’iJ  e 

(55c) 

QtjrQcj+cicje  MZ 

(55d) 

The  only  z  dependence  is  in  the  exponential  factor  so  that  far  downstream  all 
the  functions  approach  a  constant  value.  In  general  closed  form  solutions 
to  the  system  of  equations  (Eqs.  (54))  cannot  be  obtained.  There  are  several 
approximate  methods  for  solving  systems  of  linear  ordinary  differential 
equations  with  variable  coefficients  given  in  Chapter  7  of  Ref.  10.  These 
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methods  generally  result  in  series  expansions  in  powers  of  2  .  The  perturba¬ 
tion  quantities  will  have  boundary  conditions  applied  at  the  actuator  (2  =  0) 
and  at  downstream  infinity  (2->oo).  A  series  expansion  in  terms  of  powers 
of  z  is  not  particularly  useful  for  applying  the  downstream  infinity  boundary 
conditions.  However,  because  of  the  special  form  of  the  coefficients  in 
Eqs.  (54)  we  may  use  Theorem  8.1  of  Reference  11  to  construct  asymptotic  solu¬ 
tions  for  large  z  .  Essentially  the  referenced  theorem  states  that  the  solu¬ 
tion  for  a  system  of  equations  whose  coefficients  are  asymptotically  equal  to 
a  constant  (as  noted  from  Eqs.  (55))  is  asymptotically  equal  to  the  solution  of 
the  syjtem  with  constant  coefficients.  These  asymptotic  solutions  are  most 
easily  found  directly  from  Eqs.  (55).  Let  the  asymptotic  solution  for  large 
2  for  the  perturbation  amplitude  functions  be  denoted  with  a  superscript 
star.  Then  these  solutions  will  be  of  the  form 


v/  - 

(56a) 

«sT  -  C2€AZ- 

(56b) 

where  A  satisfies 

a*  »  c3ex* 

(56c) 

+  K  X) 

(?/(2x  +  J^X4) 

- 0 

74,  *♦*’,>*) 

&3Z.+  +  ^32  A  t  ?0 

(^3S->  #33  At  P33\+  A) 

(57) 

Eq.  (57)  is  obtained  by  direct  substitution  of  Eqs.  (56)  into  Eqs. 
(54).  Eq.  57  is  an  seventh  order  polynomial  in  X  with  complex  valued  coef¬ 
ficients.  After  considerable  effort  no  analytic  solutions  to  this  equation 
were  found  and  it  was  concluded  that  the  equation  would  have  to  be  solved 
numerically.  Assuming  these  roots  are  known  the  analysis  proceeds  as  follows. 

The  first  step  in  completing  the  solution  is  to  put  Eqs.  (54)  into 
standard  fi-st  order  form.  This  is  accomplished  by  differentiating  Eq.  (54a) 
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>t  II  III 

and  obtaining  simultaneous  solutions  for  V2  ,  U0  and  Uli  from  the  resulting 
three  equations.  This  results  in  the  following  system  of  equations. 

A  =  /,  Vz  4  Az  A  4  Atz  O-.  *  1,4  Uo  4  As  A  4  At  *A,i  A  +  -A, a  U,V  (58a) 

<4*  =  A,  A  4  Aiz  V'  +123  U-.  +  A4U' t  As  U  +/„<*,  Am  4  (58b) 

<  *  A,  vz  +  l3i  v2 '  + 133  +  J,4  <±' + As  4 4  -4<.a, 4  Ai  A  4  As  <4  (58c } 

The  Ay's  have  been  determined  but  they  will  not  be  given  at  the  present  time 
because  they  would  overly  complicate  this  presentation.  The  conversion  of 
Eqs.  (54)  to  standard  form  is  completed  by  renaming  the  dependent  variables 


as  follows 


A  '■  V 
4  -  ^ 

4 

A =  is 

a<  =  tffr 

ir 

• ». 


At  this  time  we  introduce  matrix  notation.  Let  (^-denote  the  column  matrix 
[  Lj..  ] .  Then  Eqs.  (58)  plus  the  definition  of  the  lead  to  the  following 
system  of  equations 

(S9) 

where  ^  ^  ]  and  A  is  an  eight  by  eight  square  matrix  given  by 
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(60) 


The  planned  method  for  determination  of  is  to  assume  that  y,  is  of  the  form 

1  =  <6l) 
where  [  D  ]  is  a  diagonal  matrix  with  the  asymptotic  solutions  for  the 
along  the  diagonal  and  |  is  a  column  matrix.  Under  these  assumptions, 

Eq.  (59)  becomes 

O.'  =  tfi]  £  (62) 


where 


:bj  =  [f0f’  {[A]  Ojil  -  Lf'z} 


Now  the  power  series  approximation  methods  given  in  Ref.  10  may  be 
applied  to  Eq.  (62)  to  solve  for  .  The  assumed  form  of  the  solution  for 
is  essentially  the  asymptotic  solution  for  large  e  times  a  series  expansion 
in  powers  of  z  .  The  rationale  for  choosing  this  representation  was  that  the 
asymptotic  solution  should  dominate  the  behavior  for  large  z  and  the  series 
solution  should  give  accurate  results  for  small  z  such  that  the  composite 
solution  should  give  reasonably  accurate  results  near  the  points  a  =  0  and 
z-»«o  where  the  boundary  conditions  are  to  be  imposed. 

The  planned  method  of  solution  became  unduly  complicated  when  we 
were  unable  to  analytically  determine  the  values  of  A.  that  satisfied  Eq.  (57). 
Numerical  solution  of  this  equation  was  judged  not  to  be  a  promising  approach 
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for  the  following  reasons.  The  coefficients  of  the  powers  of  X  appearing  in 
Eq.  (57)  are  functions  of  C  »C^+jC£  ,  the  complex  damping  factor  of  the 
assumed  disturbance  form,  which  is  unknown  at  this  stage  of  the  analysis. 
Therefore  Eq.  (57)  would  have  to  be  solved  in  conjunction  with  the  matching 
conditions  between  upstream  and  downstream  perturbation  quantities.  These 
matching  conditions  would  be  of  the  same  nature  as  those  used  in  the  two-dimen¬ 
sional  problem.  They  would  be  applied  on  a  radially  integrated  basis  across 
the  annulus  as  was  done  for  the  steady-flow  problem  in  Ref.  1  and  integral 
moments  of  the  three  fundamental  matching  conditions  would  be  used  to  generate 
a  sufficient  number  of  equations  to  make  the  problem  determinate.  However, 
since. we  don't  know,  apriori,  how  many  of  the  X.s  which  satisfy  Eq.  (57)  have 
negative  or  zero  values  for  their  real  parts  it  is  not  certain  how  many  match¬ 
ing  conditions  are  required.  This  type  of  uncertainty  in  the  formulation  can 
be  handled  but  would  require  an  extensive  effort  to  computerize  the  method 
of  solution.  At  this  point  it  was  decided  that  the  amount  of  required  numerical 
effort  would  subvert  any  physical  interpretation  of  the  results  and  the  sta¬ 
bility  analysis  was  discontinued.  It  appears  that  the  stability  of  the  three- 
dimensional  problem  would  be  best  studied  from  the  response  viewpoint. 

D.  CONCLUDING  REMARKS  ON  THEORETICAL  INVESTIGATIONS 

During  this  program,  the  two-dimensional  stability  theory  for  wholly 
subsonic  compressible  flows  was  modified  through  the  revision  of  the  energy 
matching  conditions.  Also  the  formulation  of  the  theory  to  supersonic  flow 
has  been  initiated.  A  stability  theory  for  three  dimensional  flow  was  developed 
but  not  implemented  because  of  extensive  computational  requirements.  It 
appears  that  the  stability  of  three-dimensional  flows  would  be  best  studied 
through  analysis  of  the  response  of  the  blade  row  rather  than  a  stability 
theory  approach. 
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SECTION  III 


EXPERIMENTAL  ROTATING  STALL  RESEARCH 

A.  INTRODUCTION 

As  a  part  of  the  work  under  a  previous  program  Contract  AF  33(615)- 
3357,  an  annular  cascade  facility  was  designed  and  fabricated.  Its  principal 
purpose  is  to  provide  detailed  fundamental  experimental  data  during  and  prior 
to  the  occurrence  of  rotating  stall  in  order  to  improve  our  understanding  of 
the  phenomena  and  for  use  as  a  guide  in  improving  the  theoretical  analysis. 

The  facility  has  also  been  used  to  evaluate  the  operation  of  a  prototype 
rotating  stall  control  system  (Reference  12)  and  to  provide  acoustic  data  for 
comparison  with  theory  (References  13  and  14) .  This  section  presents  the 
results  of  fundamental  experiments  on  rotating  stall. 

In  previous  programs  using  the  annular  cascade  facility  (References 
1  and  2),  rotating  stall  has  been  induced  on  a  variety  of  stator  rows  with 
differing  geometries,  and  on  a  rotor  row  with  and  without  circumferential 
inlet  distortion.  In  addition,  a  limited  examination  of  the  effects  of 
rotor-stator  interference  on  the  inception  of  rotating  stall  was  included  in 
Reference  1.  The  latter  study  showed  that  the  effect  of  close  coupling  be¬ 
tween  the  rotor  and  stator  in  a  stage  can  lead  to  significant  delay  in  rotating 
stall  inception  over  that  which  occurs  on  an  isolated  rotor.  Moreover,  the 
delay  in  rotating  stall  inception  increased  as  the  stator  stagger  angle  was 
increased,  that  is  as  the  stator  row  was  unloaded.  A  preliminary  correlation 
was  made  between  these  experimental  results  and  the  predictions  from  a  two- 
dimensional  rotating  stall  stability  theory  for  a  stage.  The  theoretical 
results  also  showed  a  delay  in  inception,  but  the  numerical  predictions  of 
the  theory  were  generally  inaccurate  due  to  a  lack  of  steady-state  performance 
data  for  the  rotor  and  stator  rows  in  the  stage.  The  theory  requires  these 
steady-state  data  as  inputs. 
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The  major  portion  of  the  current  experimental  study  was  designed  to 
provide  basic  information  on  rotor-stator  interference  and  to  provide  the 
steady-state  performance  data  required  by  the  rotating  stall  stability  theory. 
In  addition,  a  limited  study  was  performed  to  determine  the  effect  of  circum¬ 
ferential  inlet  distortion  on  rotating  stall  inception  in  the  closely  coupled 
stage .  The  results  of  these  studies  have  been  used  for  comparison  with  the 
predictions  of  the  rotating  stall  stability  theory  for  multiple  blade  rows. 
Moreover,  the  experiments  have  provided  some  unexpected  results  which  are 
important  in  their  own  right. 

All  of  the  above  rotating  stall  investigations  were  performed  in  a 
configuration  of  the  annular  cascade  with  high  hub-to-tip  ratio  (0.80). 

This  hub-to-tip  ratio  was  chosen  originally  because  the  data  would  be  used 
to  guide  and  verify  the  development  of  a  rotating  stall  theory  that  was 
basically  two-dimensional.  However,  rotating  stall  frequently  originates 
on  the  front  stages  of  axial  flow  compressors  where  the  hub-to-tip  ratio  is 
usually  quite  low  and  three-dimensional  effects  may  be  important.  Thus, 
during  the  latter  portion  of  the  current  program,  a  major  modification  was 
made  to  the  annular  cascade  facility  to  provide  it  with  a  low  hub-to-tip 
ratio  (0.44)  configuration.  The  final  tests  in  the  current  program  consisted 
of  checkout  of  the  new  configuration  and  determination  of  the  rotating  stall 
inception  boundary  on  a  low  hub-to-tip  ratio  rotor. 

The  order  of  presentation  of  the  experimental  program  in  the 
remainder  of  this  section  is  as  follows.  A  description  of  the  high  hub-to- 
tip  ratio  configuration  of  the  annular  cascade  facility  is  presented  in 
Section  III-B.  Section  III-C  presents  the  results  of  the  study  of  rotor- 
stator  interference  effects  with  and  without  circumferential  inlet  distortion. 
Section  II1-D  presents  a  description  of  the  low  hub-to-tip  ratio  configuration 
of  the  annular  cascade  and  the  preliminary  rotating  stall  inception  tests  on 
an  isolated  rotor.  Finally,  in  Section  III-E,  a  summary  of  the  experimental 
program  is  presented. 
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B. 


DESCRIPTION  OF  ANNULAR  CASCADE  FACILITY  WITH  HIGH 
HUB-TO-TIP  RATIO 


I .  General  Description 

The  annular  cascade  facility  consists  of  a  test  section  built  around 
the  outer  front  casing  of  a  J-79  jet-engine  compressor  with  a  Calspan  fabri¬ 
cated  hub.  The  facility  includes  a  bell-mouth  inlet  on  the  outer  casing  and 
a  bullet  nose  on  the  hub  to  provide  a  smooth  flow  of  air  to  the  test  section. 
Outlet  ducting  is  connected  to  an  independently  variable  source  of  suction  to 
provide  the  required  flow  through  the  annulus.  An  electrically  powered  two- 
speed  axial  flow  fan  is  used  as  the  source  of  suction.  Continuous  control  of 
the  mass  flow  is  achieved  through  the  use  of  variable  inlet  guide  vanes  to 
the  fan  and  a  variable  damper  in  the  fan  exit  flow.  Photographs  of  the  com¬ 
plete  annular  cascade  facility  and  of  the  test-section  portion  of  the  facility 
are  shown  in  Figures  2  and  3,  respectively. 

The  test  section  of  the  annular  cascade  forms  a  circular  annulus 
with  an  outer  diameter  of  29.35  inches  and  an  ii  mr  diameter  of  23.35  inches 
which  provides  a  hub-to-tip  ratio  of  0.80.  The  outer  casing  will  accept  up 
to  six  variable  stagger  angle  stator  rows.  The  hub  has  provision  for  two 
rotor  rows  at  the  third-  and  fifth-stage  rotor  locations  of  the  J-79  com¬ 
pressor.  Speed  control  on  each  rotor  hub  is  independently  variable  in  either 
direction  of  rotation.  The  overall  combination  of  variable  stagger  angle 
stators,  independently  powered  rotors,  and  independent  mass  flow  control 
provides  exceptional  versatility  to  the  complete  test  rig. 

The  configuration  of  the  annular  cascade  used  for  rotor-stator 
interference  studies  is  sketched  in  Figure  4.  The  sketch  shows  both  the 
third  and  the  fifth  stage  rotors.  However,  only  the  fifth  stage  rotor  was 
rotated  during  the  test  program.  The  third  stage  rotor  and  the  outer  skins 
on  the  hub  upstream  and  downstream  of  the  rotors,  and  between  the  rotors 
were  held  stationary.  The  fifth  stage  rotor  is  driven  by  a  hydraulic  motor 
mounted  inside  the  hub.  An  external  hydraulic  pump  system,  powered  by  a 
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50-horsepower  electric  motor,  is  used  to  provide  power  to  the  hydraulic  motor. 
With  this  system,  the  fifth  stage  rotor  speed  is  continiously  variable  between 
zero  and  approximately  1760  rpm.  This  rotor  hub  is  designed  to  accept 
46  blades.  The  method  of  fastening  the  blades  to  the  hub  allows  individual 
adjustment  of  the  stagger  angle.  In  addition  to  the  rotor  blades  on  the 
fifth  stage  hub,  a  symmetric  two-lobe  distortion  screen  was  mounted  on  the 
stationary  third  stage  rotor  hub  for  some  of  the  tests.  The  location  of  the 
distortion  screen  is  indicated  in  Figure  4. 

2 .  Rotor  and  Stator  lading  and  Distortion  Screen 

The  tests  performed  during  this  program  involved  either  an  isolated 
rotor  or  a  closely  coupled  rotor-stator  stage.  The  fifth  stage  rotor  and 
stator  blades  from  a  J-79  compressor  were  used  in  these  tests,  unmodified 
except  for  blade  length,  which  was  shortened  in  order  to  fit  into  the  three- 
inch  annulus  of  the  test  rig.  In  each  case,  the  spanwise  portion  of  the 
blades  closest  to  the  outer  compressor  casing  was  used.  The  set  of  rotor 
blades  has  been  designated  as  Rotor  Set  No.  1  and  the  set  of  stator  blades 
as  Stator  Set  No.  1,  in  keeping  with  the  nomenclature  used  in  References  1 
and  2.  The  geometric  characteristics  of  these  blades  have  been  presented 
in  References  1  and  2,  along  with  detailed  measurements  of  turning  and  Iojs 
performance  when  used  as  single  blade  rows  in  isolation.  The  extensive 
performance  data  are  not  repeated  here,  except  where  they  are  used  for 
purposes  of  comparison  with  the  current  test  results.  The  blade  geometric 
characteristics  are  repeated  below. 

Rotor  Set  No.  1  has  the  same  number  of  blades  and  is  at  the  same 
axial  location  as  the  fifth  stage  rotor  in  the  original  J-79  compressor. 

These  blades  have  a  NACA  65-Series  thickness  distribution  on  a  circular  arc 
mean  line . 


The  geometric  characteristics  are  listed  in  TABLE  1. 
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TABLE  1 


Geometric  Characteristics  of  Rotor  Set  No.  1 


Blade  Length 


2.98  inches 


Blade  Chord 


1.448  inches  (constant) 


Blade  Thickness;  Outer  Diameter 
Mid-Annulus 
Inner  Diameter 


0.072  inch 
0.086  inch 
0.113  inch 


Blade  Camber  Angle  (angle  between  tangents  to 


mean  camber 

line 

of  leading  and  trailing 

edges) ; 

Outer  Diameter 

23. S  degrees 

Mid-Annulus 

28.0  degrees 

Inner  Diameter 

32.7  degrees 

Blade  Twist 

See  Text 

Number  of  Blades 

46 

Solidity  at 

Mid-; 

4nnulus 

0.805 
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The  blades  in  Rotor  Set  No.  1,  have  a  nominal  twist  of  14.5  degrees 
over  the  three  inch  span.  However,  measurements  showed  that  in  practice,  the 
twist  over  the  center  2-1/2  inches  in  these  production  J-79  blades  varied 
between  10.8  and  14.4  degrees,  with  an  average  value  of  12.9  degrees.  The 
blades  were  installed  in  a  sequence  which  distributed  the  nonuniformities 
in  twist  over  the  complete  circumference. 

A  blade-tip  clearance  problem  was  encountered  on  installation  of 
the  rotor  assembly  in  the  annular  cascade.  The  J-79  compressor  outer  casing 
used  in  the  annular  cascade  is  about  one-tenth  inch  out  of  round  with  the 
split  line  diameter  being  the  largest.  In  the  final  assembly,  the  rotor  was 
mounted  in  the  casing  with  shims  along  the  split  line  which  provided  the 
following  tip  clearances:  0.030  inch  on  the  top  and  bottom,  and  0.046  to 
0.048  inch  on  the  sides.  Since  the  stagger  angles  of  the  rotor  blades  are 
adjustable,  the  quoted  tip  clearances  apply  only  to  the  mid-chord  pivot 
points.  Clearance  at  the  blade  leading  and  trailing  edges  vary  slightly  with 
changes  in  stagger  angle  from  the  angle  for  which  the  tip  contour  was  deter¬ 
mined.  The  stagger  angle  used  to  determine  the  tip  contour  was  the  maximum 
which  was  expected  to  be  used.  The  choice  of  lower  reference  stagger  angles 
would  have  resulted  in  interference  between  the  blades  and  outer  casing  at 
stagger  angles  much  larger  than  the  reference. 

Stator  Set  No.  1,  has  the  same  number  of  blades  and  is  located  at 
the  same  axial  location  as  the  fifth  stage  stator  row  in  the  original  J-79 
compressor.  Its  geometric  characteristics  are  listed  in  Table  2,  and  the 
axial  location  is  shown  in  Figure  4. 

For  some  of  the  tests,  a  circumferential  distortion  screen  was 
mounted  on  the  third  stage  rotor  hub,  upstream  cf  the  rotor-stator  stage. 

The  screen  has  a  symmetric  two-lobed  pattern  with  each  lobe  encompassing  a 
circumferential  extent  of  60  degrees.  The  lobes  of  the  screen  are  perforated 
steel  plates.  Detailed  calibration  of  the  distortion  generated  by  the  screen 
has  been  presented  in  Reference  1.  Briefly,  the  total  pressure  deficit  in 


TABLE  2 


Geometric  Characteristics  of  Stator  Set  No.  1 


Blade  Length 

2.98  inches 

Blade  Chord; 

Outer  Diameter 

Inner  Diameter 

1.316  inches 

1.290  inches 

Blade  Thickness; 

Outer  Diameter 

Inner  Diameter 

0.1524  inch 

0.1132  inch 

Blade  Camber  .Angle  (angle  between  tangents  to 
mean  camber  line  of  leading  and  trailing  edges) 

35.6  degrees 

Outer  Diameter 

0  degree 

Blade  Twist; 

Mid -Annulus 

0.28  degree 

positive  twist 

Inner  Diameter 

1.42  degree 

reduces  stagger 

S4 

0.85 


Number  of  Blades 
Solidity  at  Mid-Annulus 
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the  wake  of  each  lobe  was  approximately  constant  over  a  60  degree  circumfer¬ 
ential  extent.  The  magnitude  of  the  total  pressure  drop  in  the  screen  wake 
was  approximately  1.1  times  the  dynamic  pressure  measured  far  upstream  of 
the  screen. 


3 .  Instrumentation 

The  outer  casing  of  the  annular  cascade  test  section  has  been 
modified  to  allow  circumferential  and  radial  traverses  with  hot-wire  and 
total  pressure  probes  upstream  arid  downstream  of  the  fifth  stage  location. 
These  traverse  locations  are  shown  in  Figure  3.  Two  probes  can  be  accommoda¬ 
ted  simultaneously  in  either  of  these  traverse  mounts  with  adjustable  cir¬ 
cumferential  spacing  between  the  two  probes.  In  addition,  simple  radial 
traverses  can  be  made  at  various  other  locations  in  the  test  section.  The 
measuring  locations  used  in  this  program  are  indicated  in  Figure  4,  as 
Stations  0,  1,  and  2. 

In  all  of  the  experiments,  mass  flow  through  the  test  section  was 
determined  by  measuring  the  dynamic  pressure  in  the  constant-area  annulus 
upstream  of  the  test  section  (Station  0  in  Figure  4) . 

Total  pressure  surveys  were  made  through  traverses  with  a  seven- 
tube  total  pressure  rake  which  spanned  one-half  of  the  three-inch  annulus. 

A  complete  radial  survey  was  obtained  by  positioning  the  rake  at  two  radial 
locations,  one  adjacent  to  the  hub  and  one  adjacent  to  the  outer  casing. 

At  each  radial  location,  the  total  pressure  rake  was  aligned  with  the  flow 
by  using  a  Conrad  arrowhead  style  yawmeter  incorporated  at  the  mid-span 
location  of  the  rake.  At  axial  Station  1  in  Figure  4,  a  single  radial 
survey  was  made  midway  between  adjacent  stator  blades  to  determine  the  time 
averaged  total  pressure  rise  across  the  rotor.  At  Station  2,  downstream  of 
the  stators,  circumferential  as  well  as  radial  surveys  were  performed  to 
insure  that  the  stationary  stator  blade  wakes  were  properly  taken  into 
account.  In  all  tests,  the  pressures  detected  by  the  rake  were  photograph- 
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ically  recorded  from  a  multitube  inclined  monometer. 

Radial  surveys  of  the  temperature  rise  across  the  rotor  in  the 
rotor-stator  stage  were  made  with  a  matched  set  of  two  thermister  probes 
calibrated  to  measure  total  temperature  differences  as  small  as  0.01  degrees 
Fahrenheit,  in  an  environment  with  constant  upstream  reference  temperature. 

For  these  tests,  one  thermister  probe  was  located  at  mid-annulus  upstream 
of  the  rotor  to  provide  an  upstream  reference  total  temperature.  The  second 
probe  was  used  to  perform  radial  surveys  at  Station  1,  between  the  rotor 
and  stator. to  determine  the  radial  distribution  of  total  temperature  rise 
across  the  rotor.  Output  from  the  matched  thermister  probes  was  integrated 
for  50  seconds  on  an  integrating  voltmeter  and  the  average  of  seven  such 
integrated  readings  were  used  to  determine  a  single  data  point.  Thus,  each 
data  point  represented  a  time  average  for  350  seconds .  The  very  long 
averaging  times  were  used  in  an  attempt  to  compensate  for  random  fluctuations 
in  upstream  reference  temperature  in  the  annular  cascade.  The  annular  cascade 
uses  room  air  as  its  source  and  the  inlet  temperature  is  not  constant  in  time . 

In  addition  to  the  total  pressure  and  total  temperature  data 
measured  during  the  current  program,  velocity  and  swirl  angle  measurements 
presented  in  References  1  and  2,  will  be  used  to  complement  the  data  analysis 
in  the  presentation  of  the  results.  These  data  were  taken  at  Station  2 
(Figure  4)  with  a  linearized,  two-channel,  constant-temperature,  hot-wire 
anemometer  system  in  conjunction  with  a  crossed-wire  probe.  Readout  for 
the  hot-wire  system  was  made  by  means  of  a  two-channel  integrator-digital 
voltmeter  system.  Each  linearized  hot-wire  signal  was  integrated  for  60 
seconds  and  the  time  average  was  calculated  from  the  readings  on  the 
digital  voltmeters . 

Rotating  stall  was  detected  and  the  number  of  cells  and  propagation 
velocity  were  measured  by  using  outer-wall  static  pressure  taps  connected 
to  pressure  transducers.  Two  pressure  taps  at  different  circumferential 
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locations  and  at  an  axial  location  corresponding  to  the  quarter-chord  on 
the  rotor  blades  were  used. 

C.  EXPERIMENTS  ON  A  CLOSELY  COUPLED  ROTOR-STATOR  STAGE 

In  Reference  1,  a  limited  study  was  performed  to  study  the  effect 
of  rotor-stator  interference  on  the  inception  of  rotating  stall.  These 
experiments  were  conducted  in  the  configuration  of  the  annular  cascade 
shown  in  Figure  4.  The  experiments  were  designed  to  investigate  whether 
interference  between  a  rotor  and  a  stator  has  an  effect  on  the  conditions 
for  inception  of  rotating  stall.  The  rotor  row  and  the  stator  row  us  id 
in  these  experiments  had  previously  been  tested  in  isolation  so  that  the 
conditions  under  which  the  individual  blade  rows  undergo  inception  were 
known.  The  inlet  flow  conditions  relative  to  each  blade  row  were  such 
that,  in  the  absence  of  interference,  rotating  stall  would  occur  first  on 
the  rotor.  These  studies  were  performed  at  a  fixed  setting  of  the  rotor 
blade  stagger  angle,  but  with  four  different  settings  of  the  stator  blade 
stagger  angle.  The  combined  results  showed  an  increasingly  larger  delay  in 
rotating  stall  inception  in  comparison  to  the  isolated  rotor  as  the  stator 
stagger  angle  was  increased,  that  is  as  the  stator  row  was  unloaded.  Thus, 
it  was  concluded  that  rotor-stator  interference  does  have  an  effect  on 
rotating  stall  inception  and  that  the  magnitude  of  the  effect  depends  on 
the  loading  of  the  blade  row  which  would  be  unstalled  in  the  absence  of 
interference . 

An  attempt  was  made  to  correlate  the  experimental  inception  points 
on  the  stage  with  the  predictions  of  the  two-blade  row  rotating  stall 
stability  theory.  The  results  of  the  correlations  were  encouraging  in  that 
the  theory  predicted  that  blade  row  coupling  does  stabilize  the  stage. 

However,  the  numerical  predictions  were  inaccurate,  due  to  a  lack  of  steady 
state  performance  data  for  the  blade  rows.  The  theory  requires  these  steady 
state  performance  data  as  inputs.  The  major  portion  of  the  current  experimen¬ 
tal  program  was  designed  to  extend  the  experimental  data  so  that  valid  com¬ 
parisons  could  be  made  between  theory  and  experiment. 


The  steady-state  blade  row  performance  data  required  as  inputs  to 
the  rotating  stall  stability  theory  are  as  follows:  The  relative  losses 
through  the  rotor  row,  the  turning  angle  of  the  flow  through  the  rotor, 
the  losses  through  the  stator  row,  and  the  turning  angle  of  the  flow  through 
the  stator  row.  Each  of  these  quantities  will  be  addressed  in  turn  in  the 
following  discussion.  A  sketch  showing  the  velocity  notation  is  presented 
in  Figure  5 . 


1 .  Relative  Losses  Through  the  Rotor 

The  basic  experimental  data  required  as  an  input  to  the  rotating 
stall  theory  are  the  steady-state  losses  across  the  compressor  blade  rows. 
These  are  required  in  a  coordinate  system  fixed  to  the  blade  row  (relative 
losses) .  The  determination  of  the  relative  losses  through  a  rotor  row 
requires  the  measurement  of  total  pressure  rise  in  a  duct-fixed  coordinate 
system  and  a  measurement  of  the  work  done  by  the  rotor.  In  previous  tests 
on  an  isolated  rotor,  the  work  done  by  the  rotor  was  determined  from  the 
velocity  components  measured  downstream  of  the  rotor  with  a  crossed  hot-wire 
anemometer  system.  However,  in  the  current  program,  the  tests  were  performed 
on  a  closely  coupled  rotor-stator  stage  and  the  measurement  of  the  components 
of  velocity  between  the  rotor  and  stator  was  not  possible.  Thus,  an  attempt 
was  made  to  determine  the  rotor  work  by  measuring  the  total  temperature  rise 
across  the  rotor  instead  of  the  velocity  components.  It  will  be  demonstrated 
that  the  rotor  work  coefficients  calculated  from  the  total  temperature 
measurements  between  the  rotor  and  stator  are  reasonably  consistent  in  the 
mid-annulus  region  with  the  work  coefficients  calculated  from  velocity 
measurements  downstream  of  the  isolated  rotor.  However,  the  total  temperature 
results  are  not  accurate  enough  to  provide  an  estimate  of  the  work  done  by 
the  rotor  for  use  in  calculating  the  relative  losses  across  the  rotor  in  a 
rotor-stator  stage.  On  the  other  hand,  application  of  the  rotating  stall 
stability  theory  still  requires  the  relative  losses  as  an  input.  In  view  of 
the  disappointing  temperature  results,  the  rotor  work  determined  from  the 
isolated  rotor  tests  will  be  used  in  combination  with  the  total  pressure  rise 


across  the  rotor  in  the  rotor-stator  stage  to  estimate  the  relative  rotor 
loss  coefficient.  The  experimental  results  and  the  method  of  data  analysis 
are  described  below.  The  total  pressure  measurements  are  presented  first. 

Total  pressure  surveys  between  the  rotor  and  stator  (Station  1  in 
Figure  4)  were  performed  for  a  rotor  stagger  angle, <5„M,  of  40  degrees 
measured  at  mid-annulus  and  for  three  different  stator  stagger  angles,  <5S„  = 
28.2,  37.2  and  47.2  degrees,  again  measured  at  mid-annulus.  The  mean  axial 
velocity,  U0  ,  measured  far  upstream  of  the  stage  was  held  constant  at 
60  ft/sec  in  these  tests  and  measurements  were  made  for  a  wide  range  of  rotor 
speeds.  With  the  two  smaller  stator  stagger  angles,  rotating  stall  inception 
was  sudden  and  of  large  amplitude.  In  these  cases,  rotating  stall  displayed 
hysteretic  behavior  at  inception.  That  is,  as  rotor  rpm  was  increased 
rotating  stall  occurred  at  a  higher  value  of  rpm  than  the  value  at  which  it 
disappeared  on  decreasing  the  rotor  rpm.  The  stage  with  <5Sm  =  47.2  degrees 
did  not  display  hysteresis  in  rotating  stall  inception.  In  this  case,  as 
rotor  speed  was  increased,  the  outer  casing  pressure  transducers  used  to 
detect  rotating  stall  indicated  increasing  turbulence  which  became  quite 
large  just  before  inception.  The  rotating  stall  which  occurred  was  inter¬ 
mittent  in  time  and  not  much  larger  in  amplitude  than  the  turbulent  pressure 
fluctuations .  Repeated  surveys  of  total  pressure  for  this  stagger  angle 
showed  that  the  data  were  not  repeatable  near  inception.  Moreover,  the 
total  pressure  rise  across  the  rotor  began  to  decrease  with  increasing  rotor 
speed  well  before  inception  occurred.  These  results  suggest  that  steady  stall 
was  occurring  on  the  rotor  well  before  rotating  stall  inception  for  a  stator 
stagger  angle  of  47.2  degrees.  The  presence  of  steady  stall  and  non- 
repeatable  total  pressures  made  it  impractical  to  continue  tests  with  this 
atypically  large  stator  stagger  angle.  Thus,  the  remainder  of  the  tests  were 
completed  only  for  stator  stagger  angles  of  28.2  and  37.2  degrees. 


45 


Radial  distributions  of  total  pressure  coefficient  between  the 
rotor  and  stator  are  shown  in  Figures  6(a)  and  6(b),  for  stator  stagger 
angles  of  28.2  and  37.2  degrees,  respectively .  The  data  obtained  in  the 
absence  of  rotating  stall  are  shown  as  open  symbols  connected  by  solid 
lines.  Data  measured  while  rotating  stall  was  occurring  are  shown  as 
variously  dashed  lines  without  symbols.  The  total  pressure  coefficient  ACp  , 

Tt 

shown  in  Figure  6,  is  defined  as 


ACP  (r) 


(64) 


where 

=  Total  pressure  far  upstream  of  the  stage  (Station  0  in 
Figure  4) . 

*r, 

=  Total  pressure  in  the  survey  plane  between  the  rotor  and 
stator  (Station  1  in  Figure  4) . 

r 

=  Radius 

u. 

=  Mean  axial  velocity  far  upstream  (Station  0) . 

P 

=  Air  density 

and 

n 

indicates  a  time  average  in  a  fixed  frame  of  reference 
downstream  of  the  rotating  blades  (this  is  assumed  to  be 
equivalent  to  a  circumferential  average  in  a  frame  of 
reference  which  rotates  with  the  blades) . 

Some  of  the  radial  surveys  were  repeated  twice  for  6Sn=  28.2  degrees 
(Figure  6(a)).  In  these  cases,  both  sets  of  data  are  plotted.  There  is 
relatively  good  agreement  between  these  repeated  data,  in  contrast  to  the 
poor  repeatability  mentioned  previously  for  <SSM  =47.2  degrees.  Note  that 
when  rotating  stall  begins,  there  is  a  drop  in  pressure  coefficient  across 
the  complete  annulus  and  that  the  stalled  pressure  coefficients  are  nearly 
independent  of  rotor  rpm  over  the  range  tested. 

The  radial  pressure  distributions  shown  in  Figure  6  were  integrated 
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radially  (area-averaged)  to  provide  an  overall  total  pressure  rise  coefficient, 

ACP  ,  across  the  rotor.  The  results  are  shown  in  Figure  7,  along  with 
T 1 

additional  data  measured  on  the  same  rotor  in  isolation,  that  is  with  the 
stator  row  removed.  The  overall  total  pressure  rise  coefficient  is  defined 
as  follows: 


ACPt  (r)  o Lr 


(65) 


where 


=  inner  wall  (hub)  radius  of  annular  cascade 
rt  =  outer  wall  (tip)  radius  of  annular  cascade 


and  ACp  Cr*> 


is  defined  in  Equation  (64) . 


Inception  of  rotating  stall  with  increasing  rotor  speed  is  indicated 
in  Figure  7  by  a  vertical  arrow.  Note  that  the  stage  stalls  at  increasingly 
higher  values  of  rotor  speed  as  the  stator  stagger  angle  is  increased.  Also, 
the  behavior  of  the  total  pressure  rise  after  stall  is  somewhat  different  on 
the  isolated  rotor  than  on  the  stage.  On  the  isolated  rotor  the  total 
pressure  rise  after  stall  increases  fairly  rapidly  with  increasing  rotor  rpm, 
while  on  the  stage  the  total  pressure  increases  only  slightly  with  increasing 
rotor  rpm. 


The  pre-stall  total  pressure  rise  across  the  isolated  rotor  is 
significantly  lower  than  the  total  pressure  rise  across  the  rotor  in  the 
stage.  At  the  time  these  data  were  first  analyzed,  it  was  speculated  that 
some  of  the  difference  may  be  due  to  differences  in  axial  location  of  the 
total  pressure  surveys  downstream  of  the  rotor.  The  isolated  rotor  data 
shown  in  Figure  7,  were  measured  prior  to  the  current  program  (Reference  2) . 
In  that  study,  the  total  pressure  surveys  were  made  approximately  3.9  inches 
(2.73  rotor  chords)  downstream  of  the  rotor-blade  trailing  edges.  In  the 
current  study  on  the  rotor-stator  stage,  the  corresponding  measurements  were 
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made  near  the  stator-vane  leading  edges,  at  an  axial  distance  of  0.70  inches 
(0.48  rotor  chords)  from  the  rotor-blade  trailing  edges.  The  larger  rotor- 
to-probe  distance  on  the  isolated  rotor  may  have  contributed  additional 
losses  in  total  pressure  from  mixing  and  from  friction  on  the  annulus  walls. 
However,  the  differences  between  the  isolated  rotor  and  the  rotor-stator 
stage  appeared  to  be  too  large  to  be  accounted  for  solely  by  mixing  and 
additional  wall  friction.  Thus,  the  total  pressure  surveys  on  the  isolated 
rotor  were  repeated  with  the  total  pressure  rake  situated  in  its  original 
downstream  location  (3.9  inches  from  the  rotor  trailing  edge)  and  also  with 
the  rake  situated  at  the  same  location  as  used  on  the  stage  (0.7  inches 
from  the  rotor  trailing  edge) .  The  downstream  rake  location  served  as  a 
check  of  the  original  isolated  rotor  data  presented  in  Figure  7. 

The  results  from  the  new  total  pressure  surveys  with  the  downstream 
rake  location  are  compared  to  the  original  (Reference  2)  surveys  in  Figure  8. 
At  rotor  speeds  up  to  850  rpm,  the  data  agree  quite  well.  At  speeds  between 
900  and  1000  rpm,  the  current  data  show  a  small  increase  near  the  hub  and 
very  small  decreases  in  the  outer  half  of  the  annulus.  Area-weighted 
averages  of  the  new  profiles  provide  almost  identical  values  of  A  CPr  as 
those  found  from  the  Reference  2  data,  even  for  the  900,  950  and  1000  rpm 
cases.  Thus,  except  for  a  slight  radial  redistribution  of  the  profiles,  the 
current  data  agrees  quite  well  with  the  original  data.  The  radial  redistri¬ 
bution  may  have  been  caused  by  small  geometric  differences  in  the  hub  of 
the  annular  cascade,  which  was  rebuilt  between  the  two  sets  of  surveys. 

Radial  distributions  of  total  pressure  coefficient  measured  during 
the  current  program  with  the  rake  close  to,  and  far  from,  the  isolated  rotor 
are  compared  in  Figure  9.  Here,  the  data  measured  close  to  the  rotor  are 
consistently  a  small  amount  larger  than  the  data  measured  far  from  the  rotor. 
This  trend  extends  over  most  of  the  annulus  and  occurs  at  low,  as  well  as 
high  rotor  rpms.  However,  the  apparent  increase  in  total  pressure  close  to 
the  isolated  rotor  is  only  a  fraction  of  the  increase  in  total  pressure  across 
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the  same  rotor  in  a  rotor-stator  stage.  This  is  illustrated  in  Figure  10, 
where  the  radial  distributions  of  total  pressure  coefficients  for  the  rotor 
in  the  stage  prior  to  rotating  stall  inception  are  compared  to  similar  data 
for  the  isolated  rotor.  In  this  figure,  all  of  the  data  were  measured  at 
the  same  axial  location  0.7  inches  from  the  rotor  trailing  edge.  The  total 
pressure  coefficients  for  the  rotor  in  the  stage  are  larger  than  those  for 
the  isolated  rotor  and  the  differences  increase  with  rotor  speed  up  to 
1000  rpm.  At  1100  rpm,  only  the  stage  with  <5SM  =  37.2  degrees  remained 
unstalled,  so  a  comparison  is  not  possible  for  this  speed. 


All  of  the  total  pressure  data  presented  thus  far  are  summarized 
in  Figure  11,  where  area-averaged  total  pressure  rise  coefficients  based  on 
relative  inflow  velocity  to  the  rotor  are  shown  as  a  function  of  flow  co¬ 
efficient.  The  notation  in  this  figure  is  as  follows: 


Pressure  Coefficient,  AC, 


Work  Coefficient, 


(Pr.-PrJ  _  (Pr,  -  PT.) 


2  W, 


% 


U. 


(66) 


(67) 


Flow  Coefficient, 
where : 


C?  =  -ar- 
Wb 


,  p  g  . 

Cj.  =  Dynamic  pressure  corresponding  to  a  mean  axial  velocity,  ( y  Ut  j 

cf.i  =  Dynamic  pressure  corresponding  to  mean  velocity  relative 
to  the  rotor,  -j-  ( uf  +  wj  ) 
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P  =  Time-averaged  total  pressure  at  mid-annulus  upstream  of  rotor 
(Measured  in  absolute  coordinate  system) 

PT  =  Time-averaged  total  pressure  at  given  radius  downstream  of  rotor 
(Measured  in  absolute  coordinate  system) 

(Pr>-  PT)  -  Area  weighted  average  of  total  pressure  rise  across  rotor 

U0  =  Mean  axial  velocity  at  mid-annulus  far  upstream  of  rotor 

W,  =  Circumferential  flow  velocity 

(Measured  in  absolute  coordinate  system  downstream  of  isolated 
rotor) 

u/(  =  Time  average  of  W, 

wb  =  Rotor  Velocity  (Constant  in  time) 

Wb  =  Area  weighted  average  of  rotor  velocity 
WlW,  =  Area  weighted  average  of  Wb  W, 

The  area  weighted  averages  were  computed  by  numerically  integrating 
the  various  quantities  over  the  span  of  the  annular  cascade  in  the  same 
fashion  as  indicated  by  Equation  65 . 

As  can  be  seen  in  Figure  11,  the  current  and  original  (Reference  2) 
data  measured  with  the  total  pressure  rake  3.9  inches  downstream  of  the 
isolated  rotor  agree  very  well.  The  current  data  measured  with  the  rake  0.7 
inches  from  the  isolated  rotor  fall  a  small,  almost  constant,  amount  above 
the  downstream  measurements.  The  difference  in  this  case  may  well  be  due 
to  increased  friction  and  mixing  losses.  However,  there  is  a  significant 
increase  in  total  pressure  rise  across  the  rotor  when  the  stator  row  is 
added  immediately  downstream  of  the  rotor.  This  is  true,  even  when  the 
total  pressure  rake  is  situated  in  the  same  location  for  measurements  on  the 
isolated  rotor  as  it  is  for  measurements  on  the  stage.  Thus,  the  presence 
of  the  stator  row  appears  to  have  a  beneficial  upstream  influence  on  the 
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performance  of  the  rotor.  It  is  not  known  whether  the  presence  of  the  stators 
increases  the  work  done  by  the  rotor  or  decreases  the  losses  through  the  rotor. 
A  work  coefficient  calculated  from  hot-wire  measurements  downstream  of  the 
isolated  rotor  is  shown  in  Figure  11,  for  comparison.  If  this  work  co¬ 
efficient  applies  to  the  staged  rotor  as  well  as  the  isolated  rotor,  then 
the  differences  between  this  curve  and  the  total  pressure  rise  data  give  the 
relative  loss  coefficient  across  the  rotor.  As  noted  at  the  beginning  of 
this  section,  one  of  the  inputs  required  by  the  rotating  stall  stability 
theory  is  the  relative  loss  coefficient  across  the  rotor.  Thus,  an  attempt 
was  made  to  determine  the  actual  work  coefficient  of  the  rotor  in  the  stage 
by  measuring  the  total  temperature  rise  across  the  rotor.  The  results  are 
described  below. 

The  relative  loss  coefficient,  (-A Co  ),  across  the  rotor  is 

rT  a 

given  by  the  difference  between  the  work  done  by  the  rotor  and  the  total 
pressure  rise  generated  by  the  rotor.  For  incompressible  flow,  this  is 


(~ACP  ) 

•  o 


2  Cp  A  Tr  _ 

=  - -  AC, 


(68) 


where  Cp  =  specific  heat  at  constant  pressure  (Cp  =  6006  -f tVscc"#  for  air) 


AC 


pr  =  time  averaged  total  pressure  rise  coefficient  across  rotor 


<Pt,  -  t 


in  absolute  coordinate  system 


Pr  -  time  averaged  total  pressure  (pounds/ft  )  in  absolute 
coordinate  system 

ATt  =  time  averaged  total  temperature  rise  (  fT  -Tr  )  across  rotor 
1  in  absolute  coordinate  system  ’  * 

Tt  =  time  averaged  total  temperature  (deg .  Rankine)  in 
absolute  coordinate  system. 

U„  =  mean  axial  velocity  (ft/sec)  far  ahead  of  rotor 

p  =  density  of  air  (slugs/ft^) . 
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and  subscripts 


0  -  denotes  conditions  far  upstream  of  rotor 

1  -  denotes  conditions  downstream  of  rotor,  but 
ahead  of  stator  vanes 

R  -  denotes  conditions  in  coordinate  system  fixed 
to  rotor,  (relative  coordinate  system) . 


The  relative  loss  coefficient  can  also  be  expressed  in  terms  of  the  swirl 
velocity  downstream  of  a  rotor.  In  this  case 


2  VJb  (W,  - 


-  AC, 


(69) 


In  these  tests,  VJ  =  O 


so  that 


(-AC»  )  = 


2  W, 


-  AC. 


The  total  temperature  rise,  atT)  ,  measurements  were  made  with 
a  matched  set  of  two  thermister  probes  as  described  in  Part  A.  3  of  this 
section.  Radial  distributions  of  the  rotor  work  coefficient  T,i. 

u. 

calculated  from  some  of  the  total  temperature  rise  surveys  are  compared  in 


section.  Radial  distributions  of  the  rotor  work  coefficient  - -  f  7 

calculated  from  some  of  the  total  temperature  rise  surveys  are  compared  in 
Figure  12,  to  the  total  pressure  rise  coefficient,  ACp^  ,  and  to  the 
rotor  work  coefficient,  ^  ,  measured  with  the  stator  row  removed.  The 

data  in  Figure  12(a)  are  for  a  stator  stagger  angle  of  28.2  degrees.  Results 
from  temperature  surveys  performed  on  two  different  dates  are  shown  for  two 
rotor  speeds.  As  can  be  seen,  there  is  considerable  scatter  in  the  tempera¬ 
ture  survey  results  repeated  on  different  dates.  Moreover,  at  the  lower 
rotor  speed,  (700  RPM)  the  rotor  work  coefficient  is  smaller  than  the  total 
pressure  rise  coefficient  over  most  of  the  annulus.  This  would  indicate  a 
negative  relative  loss  coefficient,  which  is  not  possible.  Results  are 
shown  for  a  stator  stagger  angle  of  37.2  degrees  in  Figure  12  (b) .  In  this 
case,  the  temperature  measurements  were  not  repeated  on  different  dates. 
However,  the  trends  in  the  data  are  similar. 


The  very  large  rise  in  the  apparent  work  coefficient  near  the  hub, 
^r»0,  is  believed  to  be  a  result  of  heat  added  to  the  flow  by  the 
hydraulic  motor  used  to  drive  the  rotor.  The  latter  is  situated  inside  the 
hub  and  becomes  hotter  as  the  test  progresses.  The  addition  of  external  heat 
from  the  motor  would  explain  a  peculiarity  of  the  temperature  surveys  noted 
during  the  tests.  It  was  noted  that  measurements  near  the  hub  showed  a  sub¬ 
stantially  smaller  temperature  rise  at  the  beginning  of  a  test  than  a 
similar  measurement  made  near  the  end  of  the  test.  On  the  other  hand, 
repeated  temperature  tests  at  mid-annulus  did  not  exhibit  this  consistent 
behavior.  Moreover,  temperature  surveys  at  very  low  rotor  and  flow  speeds 
also  showed  a  temperature  rise  near  the  hub  when  it  should  not  have  been 
possible  to  detect  such  a  rise.  Thus,  it  was  concluded  that  the  temperature 
rise  data  measured  near  the  hub  were  contaminated  by  the  addition  of  external 
heat . 

In  summary,  the  above  comments  indicate  that  the  total  temperature 
results  are  not  accurate  enough  to  provide  an  estimate  of  the  relative  losses 
across  the  rotor  in  a  rotor-stator  stage.  On  the  other  hand,  application  of 
the  rotating  stall  stability  theory  still  requires  the  relative  losses  as  an 
input.  In  view  of  the  disappointing  temperature  results,  the  rotor  work  co¬ 
efficient,  ,  determined  from  the  isolated  rotor  tests  was  used  in 

combination  with  the  total  pressure  rise  across  the  rotor  in  the  rotor-stator 
stage  to  estimate  the  relative  rotor  loss  coefficient,  (see  Equation  70) . 

The  radial  distributions  of  this  work  coefficient  are  shown  in  Figures  12(a) 
and  12(b),  along  with  the  measured  total  pressure  rise  across  the  staged 
rotor,  and  the  work  coefficient  calculated  from  the  temperature  measurements. 
It  is  worth  noting  that  even  though  the  temperature -derived  work  coefficients 
are  not  very  accurate,  they  do  appear  to  be  reasonably  consistent  in  the  mid¬ 
annulus  region  with  the  velocity-derived  isolated  rotor  work  coefficients. 

The  data  in  Figure  12  show  that  at  low  rotor  speeds,  the  isolated 
rotor  work  coefficient  is  nearly  identical  to  the  total  pressure  coefficient 
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in  the  middle  portion  of  the  annulus.  The  fact  that  the  isolated  rotor  work 
coefficient  is  very  slightly  smaller  than  the  total  pressure  rise  in  the  mid¬ 
annulus  region  should  not  be  taken  as  an  indication  of  negative  relative 
losses.  There  are  two  reasons  for  this  statement.  First,  the  hot-wire 
velocity  measurements  are  accurate  to  within  approximately  3  percent  of 
the  true  value  and  some  of  the  apparent  (small)  negative  relative  losses 
could  be  due  to  the  accuracy  of  the  hot-wire  measurements.  Second,  the  hot¬ 
wire  velocity  measurements  on  the  isolated  rotor  were  made  farther  downstream 
than  the  total  pressure  measurements  on  the  staged  rotor.  Thus,  radial  re¬ 
distribution  of  the  flow  could  account  for  the  apparent  negative  local 
relative  loss  coefficients.  In  general,  comparisons  of  the  complete  set  of 
measurements  for  all  rotor  speeds  on  the  isolated  and  staged  rotor  are 
encouraging  in  that  the  trends  shown  are  what  one  would  expect.  Even  the 
isolated  rotor  work  coefficient  after  rotating  stall  inception  appears  to 
be  consistent  with  the  data  obtained  prior  to  stall  inception.  This  is 
illustrated  in  Figure  13,  where  radial  distributions  of  work  coefficient  are 
shown  for  the  isolated  rotor  at  various  rotor  speeds,  both  unstalled  and 
stalled. 


As  can  be  seen,  the  isolated  rotor  work  coefficients  with  and  with¬ 
out  rotating  stall  appear  to  be  consistent  with  each  other,  except  for  two 
unstalled  data  points  at  1000  rpm  on  the  rotor.  For  this  rotor  speed,  radial 
distributions  of  the  isolated  rotor  work  coefficient  are  nearly  identical 
over  the  annulus  in  the  range  O  <  &<r  £  2.2  inches.  The  difference  in 
work  coefficients  for  the  stalled  and  unstalled  rotor  at  1000  rpm  occurs 
near  the  tip  where  the  unstalled  rotor  work  coefficient  decreased  signifi¬ 
cantly;  a  result  which  appears  inconsistent  with  the  remainder  of  the  un¬ 
stalled  data.  At  the  same  time,  the  stalled  rotor  work  coefcients  at  1000  rpm 
appear  consistent  with  the  unstalled  data  at  lower  rotor  speeds  and,  also 
consistent  with  the  stalled  data  at  higher  rotor  speeds.  The  reason  for  the 
inconsistent  behavior  of  the  unstalled  isolated  rotor  at  1000  rpm  is  not 
known.  However,  the  main  point  to  be  made  here  is  that  the  isolated  rotor 
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work  coefficients,  with  the  exception  of  the  unstalled  data  near  the  tip 
at  1000  rpm,  are  well-behaved  through  and  beyond  the  stall  boundary. 

Moreover,  since  they  agree  reasonably  well  with  the  work  obtained  from  the 
total  temperature  measurements  taken  on  the  rotor-stator  stage  in  the  mid¬ 
annulus  region,  the  isolated  rotor  work  coefficients  were  assumed  to  hold 
for  the  rotor-stator  stage.  Thus,  the  solid  line  in  Figure  11  was  used  as 
an  estimate  of  the  area-weighted  average  of  the  rotor  work  coefficient  for 
the  rotor-stator  stage.  The  work  coefficient  in  this  figure  was  derived 
from  the  unstalled  measurements  in  Figure  13  for  rotor  rpm's  £950(df>o.5i i) 
and  from  stalled  measurements  for  rotor  rpm's  >  950  {<$  <■  0.5n )  . 

The  area-weighted  average  of  the  rotor  loss  coefficient  was  calcu¬ 
lated  from  the  data  shown  in  Figure  11  as  follows. 
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The  results  are  shown  in  Figure  14  for  the  isolated  rotor  and  for  the  rotor 
in  the  rotor-stator  stage  with  two  different  stator  stagger  angles.  The 
data  for  the  rotor-stator  stage  at  1000  rpm  and  above,  are  based  on  the 
isolated  rotor  work  coefficients  measured  in  the  presence  of  rotating  stall. 
In  all  cases,  the  rotor  pressure  rise  coefficients,  ACP  ,  used  in  the 
calculations  were  measured  with  rotating  stall  absent.  Experimental  in¬ 
ception  of  rotating  stall  is  indicated  in  Figure  14  by  a  vertical  arrow. 

The  solid  lines  drawn  through  the  loss  data  are  spline-fitted  to  the  experi¬ 
mental  points.  These  spline-fitted  curves  were  used  as  inputs  to  the  two- 
blade-row  stability  theory  for  calculating  rotating  stall  inception. 
Application  of  the  theory  to  the  rotor-stator  stage  required  a  small  extra¬ 
polation  of  the  loss  data.  The  extrapolations  are  shown  as  dashed  extensions 
to  the  solid  lines. 


The  rotating  stall  stability  theory  requires  several  other  inputs, 
as  well  as  the  relative  losses  through  the  rotor  in  the  rotor-stator  stage. 
These  inputs  are:  the  turning  angle  of  the  flow  through  the  rotor,  the 
losses  through  the  stator  row,  and  the  turning  angle  of  the  flow  through  the 
stator  row.  The  methods  used  to  estimate  these  parameters  from  the  experi¬ 
mental  data  are  discussed  below. 


2 .  Turning  Angle  of  Flow  Through  the  Rotor 

In  previous  tests  on  an  isolated  rotor,  the  turning  of  the  flow 
through  the  rotor  was  obtained  by  measuring  the  velocity  components  down¬ 
stream  of  the  rotor  with  a  crossed  hot-wire  anemometer  system.  As  noted  at 
the  beginning  of  this  section,  it  was  not  possible  to  obtain  similar  velocity 
data  between  the  rotor  and  stator  in  the  closely-coupled  rotor-stator  stage. 
Moreover,  the  isolated  rotor  turning  data  cannot  be  used  directly,  because 
the  isolated  rotor  stalled  at  a  lower  rotor  speed  than  the  rotor  in  the 
rotor-stator  stage.  This  leaves  one  with  the  requirement  to  estimate  the 
turning  through  the  unstalled  rotor,  under  conditions  for  which  the  isolated 
rotor  was  in  rotating  stall.  In  the  previous  paragraphs,  the  work  coefficient 
measured  on  the  isolated  rotor  was  shown  to  be  insensitive  to  the  presence 
of  rotating  stall,  and  to  provide  what  appears  to  be  a  reasonable  basis  for 
calculating  the  relative  losses  through  the  rotor  in  a  rotor-stator  stage. 
Thus,  the  isolated  rotor  work  coefficient  was  also  used  to  estimate  the 
effective  turning  of  the  flow  through  the  rotor. 


For  incompressible  flow,  the  time-averaged  absolute  outlet  swirl, 
tan  /3,  ,  from  the  rotor  is  given  by  (see  Figure  5) , 

_  W. 


U. 


and  the  time-averaged  relative  outlet  swirl,  tan/3,  ,  is  given  by 

r—.  *•  -  ” 
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Note  that  in  the  experiments,  Wfc  was  held  constant  in  time,  so  Wk  «  . 


Area-weighted  averages 


of  the  absolute  and  relative  swirls  are  given  by 


and 
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From  continuity,  the  area-weighted  axial  flow  velocity,  0,  ,  is  equal  to 
the  mean  axial  flow  velocity,  U,  ,  far  upstream  of  the  rotor.  It  was 
assumed  that 


and 
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With  these  assumptions,  one  has 
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The  isolated  rotor  work  coefficient,  (2  /  U0*  )  ,  has  already  been  cal¬ 

culated  as  a  function  of  rotor  rpm.  (Figure  11  presents  a  slightly  different 
form  of  this  coefficient  as  a  function  of  flow  coefficient.)  The  area 
average  Wb/  U„  is  easily  calculated  from  the  annulus  geometry  and  rotor  rpm. 

The  results  from  Equations  72  and  73  are  plotted  in  Figure  15,  as 
a  function  of  rotor  rpm.  Note  that  the  negative  of  the  relative  swirl (- to**  /3,^) 
is  plotted  to  present  both  curves  on  the  same  graph.  Also  shown  in  Figure  15, 
are  the  results  from  direct  measurements  of  tan  and  -tan  /3t#  which  were 
reported  in  Reference  2.  In  the  absence  of  rotating  stall  (rpm  <  1000), 
the  swirl  calculated  from  the  rotor  work  agrees  very  well  with  the  direct 
measurements.  In  the  stalled  region  (rpm  >  1000),  the  direct  measurements 
of  swirl  are  less  than  the  values  calculated  from  the  rotor  work.  The  direct 
measurements  display  a  discontinuity  at  stall  inception,  whereas  the  calcu¬ 
lated  values  progress  smoothly  through  stall  inception.  Thus,  the  calculated 
swirl  values  for  the  isolated  rotor  were  used  in  the  theory  as  an  estimate 
of  the  unstalled  turning  performance  of  the  rotor  in  the  rotor-stator  stage. 


Losses  Through  the  Stator  Row 


The  losses  through  the  stator  row  in  the  rotor-stator  stage  were 
obtained  by  combining  the  results  from  total  pressure  surveys  downstream 
of  the  stators  with  the  results  from  the  total  pressure  surveys  between 
the  rotor  and  stator.  The  area-averaged  loss  coefficient  through  the  stator 
row  is  given  by 


where  ACPt^  is  the  area-averaged  total  pressure  rise  coefficient  across 
the  rotor,  (see  Figure  7)  and  A^r^  is  the  area-averaged  total  pressure 
rise  coefficient  across  the  stage. 


58 


The  area-averages  of  the  total  pressure  rise  coefficients  <iCp_ 
and  ACp„  were  obtained  in  a  different  manner  for  each  axial  location.  At 

Ti 

axial  location  1,  between  the  rotor  and  stators,  a  single  radial  survey  was 
made  mid-way  between  adjacent  stator  blades  to  determine  the  time-averaged 
total  pressure  rise  across  the  rotor.  These  time-averaged  data  were  then 
integrated  radially  to  obtain  an  area  average  (Equation  65).  At  location  2, 
downstream  of  the  stators,  it  was  necessary  to  perform  circumferential  as 
well  as  radial  surveys  to  ensure  that  the  stationary  stator  blade  wakes  were 
properly  taken  into  account.  In  this  case,  radial  surveys  were  performed 
at  7  equally-spaced  circumferential  locations  and  the  results  were  averaged 
circumferentially.  The  circumferential  averages  were  then  integrated  radially 
to  obtain  the  required  area  average.  The  circumferential  extent  of  the  surveys 
was  large  enough  to  cover  two  stator-blade  spacings .  Fourteen  radial  locations 
were  used  in  the  total  pressure  surveys  at  both  axial  locations.  Thus,  the 
area-averaged  total  pressure  coefficient  downstream  of  the  stator  represents 
an  integration  over  98  separate  spatial  locations. 


The  results  from  the  total  pressure  surveys  are  summarized  in 
Figure  16,  where  the  area-averaged  total  pressure  loss  coefficients  through 
the  stator  row  (Equation  74)  are  shown  as  a  function  of  rotor  speed.  The 
lines  drawn  through  the  experimental  data  were  spline-fitted  to  the  unstalled 
experimental  points  and  used  as  inputs  to  the  theory  to  represent  the  un¬ 
stalled  losses  through  the  stator  row. 

4 .  Turning  Angle  of  the  Flow  Through  the  Stator  Row 

The  final  input  required  by  the  rotating -stall  stability  theory  is 
the  turning  angle  of  the  flow  through  the  stator  tow.  The  prediction  of 
rotating  stall  inception  is  only  weakly  dependent  on  this  parameter,  while 
measurement  of  the  turning  through  the  stator  row  would  require  extensive 
spatial  surveys  downstream  of  the  stage  with  hot-wire  equipment.  These 
extensive  surveys  were  judge  to  be  not  necessary  in  view  of  the  weak 
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dependence  of  the  theoretical  predictions  on  the  results.  Thus,  the  turning 
through  the  stators  was  estimated  on  the  basis  of  correlations  developed  in 
Reference  1 .  The  turning  data  used  in  the  theory  has  been  presented  in 
Figure  32,  of  Reference  1  and  is  reproduced  here  in  Figure  17.  A  full  dis¬ 
cussion  of  the  correlation  is  presented  in  Reference  1. 

S .  Effect  of  Inlet  Distortion 

In  Reference  1,  the  results  of  an  extensive  series  of  tests  on  an 
isolated  rotor  with  stationary  and  rotating  inlet  distortion  were  presented. 
The  tests  on  stationary  distortion  effects  on  a  rotor  row  included  both  a 
two-lobed  and  a  four-lobed  circumferential  distortion  pattern.  Although  the 
total  pressure  distortions  were  quite  high  for  both  patterns,  they  did  not 
significantly  affect  the  rotating  stall  boundary  or  the  rotating  stall 
properties  after  inception.  Similarly,  measurements  of  averaged  steady-state 
blade  row  loss  and  turning  data  prior  to  inception  showed  little  influence 
of  the  distortion.  The  tests  with  moving  distortion  were  performed  by  ro¬ 
tating  the  two-lobed  pattern  about  the  compressor  axis,  both  in  the  direction 
of  blade  row  rotation  and  opposite  to  it.  It  was  found  that  the  rotating 
distortion  pattern  either  had  little  effect  on,  or  else  destructively  inter¬ 
fered  with  rotating  stall  formation,  except  when  the  pattern  was  rotated  at 
speeds  near  the  natural  stall  propagation  velocity  for  undistorted  flow. 

Under  these  latter  conditions,  the  rotor  incurred  rotating  stall  at  much 
higher  flow  coefficients  than  for  the  undistorted  flow  case.  A  linearized 
small  disturbance  analysis  was  developed  to  predict  the  effect  of  rotating 
distortion.  The  theoretical  analysis  was  in  general  agreement  with  the 
experimental  findings.  The  theory  concluded  that  if  the  steady-state  blade 
row  performance  is  unchanged  in  the  presence  of  distortion,  then  the  rotating 
stall  boundary  will  not  be  influenced  unless  the  rotating  distortion  pattern 
contains  a  Fourier  component  which  corresponds  to  the  natural  rotating  stall 
mode  for  the  system. 
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Moreover,  it  was  pointed  out  that  the  linearized  small-disturbance 
theory  would  result  in  a  similar  conclusion  for  multiple  blade  rows  such  as 
in  a  rotor-stator  stage.  That  is,  the  distortion  would  not  affect  rotating 
stall  inception  unless  the  averaged  steady-state  blade  row  performance 
changed  or  unless  the  distortion  was  rotated  at  a  speed  near  that  of  naturally 
occurring  rotating  stall.  If  the  distortion  is  stationary,  then  only  a  change 
in  averaged  steady-state  turning  and  loss  performance  would  affect  inception 
and  the  two-blade  row  rotating  stall  theory  for  undistorted  inlet  flow  should 
be  applicable  to  the  stage  with  distortion.  A  limited  investigation  of  this 
possibility  was  performed  during  the  current  program. 

The  two-lobed  circumferential  distortion  screen  used  for  tests  on 
the  isolated  rotor  in  Reference  1  was  installed  upstream  of  the  stage  as 
shown  in  Figure  4.  Detailed  calibrations  of  the  distortion  generated  by  this 
screen  are  presented  in  Reference  1.  A  brief  summary  of  the  results  is  given 
in  Part  A.  3  of  this  section. 

The  rotating  stall  inception  boundary  was  determined  on  the  stage 
with  inlet  distortion  for  stator  stagger  angles  of  28.2  and  37.2  degrees. 

The  results  are  presented  in  Table  3,  along  with  similar  results  obtained 
on  the  stage  without  distortion  and  on  the  isolated  rotor.  As  mentioned 
previously,  the  presence  of  inlet  distortion  had  very  little  effect  on 
rotating  stall  inception  on  the  isolated  rotor.  With  the  stage,  inlet 
distortion  delayed  inception  by  a  small  amount  for  both  stator  stagger  angles. 
This  unexpected  result  may  be  due  to  destructive  interference  between 
naturally  occurring  rotating  stall  on  the  stage,  and  the  disturbances  gener¬ 
ated  by  the  distortion  screen.  A  similar  result  was  observed,  in  some 
instances,  on  the  isolated  rotor  with  rotating  distortion  (Reference  1)  and 
on  a  stator  row  with  a  disturbance  generator  mounted  between  the  stator 
blades  (Reference  2) .  It  is  also  possible  that  the  delay  in  inception  was 
caused  by  changes  in  blade  row  steady-state  performance  in  the  presence  of 
inlet  distortion.  The  two  blade  row  stability  theory  has  been  used  to 
investigate  the  latter  effect.  The  experimental  steady-state  performance  data 
used  as  inputs  to  the  theory  for  this  application  are  described  below. 


61 


In  Reference  1,  it  was  found  that  the  circumferential  distortion 
screen  had  a  negligible  effect  on  the  total  pressure  rise  across  the  isolated 
rotor.  This  is  illustrated  in  Figure  18,  where  the  dimensionless  total 
pressure  rise  across  the  rotor  with  and  without  distortion  is  shown  as  a 
function  of  flow  coefficient.  The  same  result  was  assumed  to  hold  for  the 
rotor  in  the  rotor-stator  stage.  Thus,  the  undistorted  total  pressure  rises 
for  the  rotor  in  the  stage  with  the  two  stator  stagger  angles  (Figure  11) , 
were  used  to  calculate  the  corresponding  relative  loss  coefficients  for  the 
rotor  in  the  stage  with  inlet  distortion 

Also  in  Reference  1,  it  was  found  that  the  presence  of  inlet  distor¬ 
tion  had  only  a  small,  almost  constant  effect  on  the  turning  of  the  flow 
through  the  isolated  rotor.  The  results  from  that  study  are  shown  in  Figure  19. 
Here,  radial  distributions  of  average  swirl  angle  downstream  of  the  isolated 
rotor  with  inlet  distortion  are  compared  to  data  obtained  without  inlet  dis¬ 
tortion.  The  data  obtained  with  inlet  distortion  are  averages  of  20  measure¬ 
ments  made  at  different  circumferential  locations  downstream  of  the  isolated 
rotor.  In  Figure  19,  one  degree  has  been  added  to  the  data  measured  with  in¬ 
let  distortion,  so  that  direct  comparisons  can  be  made  with  the  shapes  of  the 
radial  distributions  for  the  isolated  rotor  without  distortion.  A  similar 
presentation  in  Reference  1  (Figure  20)  used  a  two  degree  increment  for  the 
distorted  flow  measurements.  However,  a  one  degree  increment  appears  to  pro¬ 
vide  a  better  correspondence  with  the  undistorted  data  at  the  higher  rotor 
speeds  near  rotating  stall  inception.  It  is  evident  from  Figure  19  that  inlet 
distortion  reduced  the  outlet  swirl  angle  from  the  rotor  by  an  approximately 
constant  one  degree  at  rotor  speeds  near  inception. 

It  was  assumed  that  the  rotor  in  the  stage  was  similarly  affected 
by  the  inlet  distortion.  This  assumption  has  two  effects  in  the  application 
of  the  theory.  The  most  obvious  effect  is  that  the  angle  of  attack  on  the 
stators  is  reduced  by  one  degree  at  any  given  rotor  speed.  The  second  effect 
is  that  the  rotor  work  coefficient  is  decreased  slightly  for  a  given  rotor 
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The  area-weighted  relative  loss  coefficient  for  the  rotor  in  the  stage  with 
distortion  then  becomes. 


9 rKOSO  KOTOR 
WITHOUT  OISTOKTIOfJ 


(76) 


where  the  last  term  on  the  right  is  a  consequence  of  the  previous  discussion 
of  Figure  18. 


In  summary  of  the  above  discussion,  the  inputs  required  for  appli¬ 
cation  of  the  stability  theory  were  obtained  as  follows.  Equations  (75) 
and  (76)  were  used  to  estimate  the  relative  losses  through  the  rotor  in  the 
rotor-stator  stage  with  distortion.  The  turning  performance  of  the  staged 
rotor  was  estimated  from  the  data  in  Figure  15,  corrected  for  one  degree 
less  turning.  Stator  row  losses  were  taken  from  the  data  in  Figure  16, 


64 


corrected  for  the  one  degree  decrease  in  angle  of  attack.  Finally,  the 
stator  row  turning  was  taken  from  the  correlation  shown  in  Figure  17. 

D.  LOW  HUB-TO-TIP  RATIO  ANNULAR  CASCADE 

At  the  conclusion  of  the  test  program  on  the  high  hub-to-tip  ratio 
configuration  of  the  annular  cascade,  major  modifications  were  made  to  the 
hub  of  the  annular  cascade  to  provide  it  with  a  lower  hub-to-tip  ratio.  The 
modified  installation  is  shown  in  Figure  20. 

The  new  hub  was  designed  to  accept  the  first  stage  rotor  from  a 
J-79  compressor  at  its  normal  and  location  relative  to  the  J-79  outer 
casing  which  forms  the  test  section  of  the  annular  cascade.  This  allows 
tests  with  the  rotor  alone,  or  with  a  stage  consisting  of  the  J-79  first 
stage  rotor  and  first  stage  stator.  The  hub  provides  a  constant  area  annulus 
with  a  hub-to-tip  ratio  of  0.44.  The  hub  has  an  elliptical  nose  fairing  at 
the  front  of  the  annulus,  and  the  outer  casing  has  the  same  bell  mouth 
inlet  as  was  used  for  the  high  hub-to-tip  ratio  configuration  (Figure  2) . 

The  rotor  uses  the  original  rotor  hub  section  from  the  J-79  com¬ 
pressor  to  retain  the  blades.  The  original  tapered  hub  was  changed  to  a 
constant  area  design  by  filling  the  inter-blade  spaces  on  the  hub  of  the 
rotor  blades.  The  tapered-end  stator  blades  were  converted  to  constant  area 
by  trimming  the  hub  ends  to  the  desired  shape.  The  configuration  of  the 
stator  blade  hub-end  fittings  was  retained  to  allow  use  of  the  original 
stator  hub  support  ring  from  the  J-79  compressor.  The  final  hub-to-tip  ratio 
(0.44)  of  the  new  configuration  was  governed  by  the  radial  length  of  the 
original  stator  blade  trailing  edges.  With  this  design,  the  rotor  blade 
stagger  angle  is  fixed,  but  the  stator  stagger  is  variable  through  the 
variable  geometry  mechanisms  on  the  J-79  compressor  casing. 
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The  new  hub  configuration  was  designed  to  incorporate  as  much  of 
the  existing  hydraulic  rotor-drive  system  and  downstream  ducting  as  possible. 
Incorporation  of  these  components  into  the  design  required  the  use  of  a  tran¬ 
sition  section  at  the  downstream  end.  The  transition,  shown  on  the  left-hand 
side  of  Figure  20,  decreased  the  annulus  area  to  the  original  high  hub-to-tip 
ratio  value.  Mass  flow  through  the  test  section  is  controlled  by  the  existing 
two-speed  axial  flow  fan  and  variable  damper  system. 

In  addition  to  the  hub  modifications  to  the  annular  cascade,  a  new 
radial  traverse  system  for  total  pressure  and  hot-wire  probes  was  built. 

The  new  traverse  is  powered  by  a  variable  speed  electric  motor  and  is  capable 
of  traversing  the  larger  radial  distances  associated  with  the  new  hub.  Micro¬ 
switches  prevent  accidental  encounters  between  the  probes  and  the  test  section 
hub  and  tip.  The  multiple-tube  total  pressure  rake  used  in  the  current  study 
has  been  modified  for  use  with  the  new  traverse.  This  rake  covers  only  a 
fraction  (1.5  inches)  of  the  total  radial  span  of  the  modified  annular  cascade. 
Full  span  surveys  are  obtained  by  traversing  the  rake  to  different  radial  loca¬ 
tions.  This  procedure  has  been  adopted  to  ensure  that  the  rake  is  aligned 
approximately  with  the  local  flow  direction,  which  is  determined  at  each  radial 
location  from  a  Conrad  arrowhead  style  yawmeter  incoiporated  on  the  rake. 

The  full  rotor-stator  stage  configuration  was  assembled  and  found  to 
fit  properly  without  interference.  Following  this,  the  stators  were  removed 
and  checkout  tests  of  the  rotor  alone  were  performed.  The  hydraulic  drive 
system  was  found  to  be  capable  of  driving  the  rotor  to  speeds  in  excess  of 
1200  rpm  without  vibration  or  interference  problems.  However,  maintaining  a 
given  mean  axial  velocity  in  the  test  section  requires  considerably  more  power 
from  the  downstream  fan  than  was  required  for  the  high  hub-to-tip  ratio  con- 
■  ration.  It  is  now  necessary  to  use  the  downstream  fan  on  high  speed  to 
•  .  the  60  feet  per  second  mean  axial  velocity  normally  used  for  the  ex- 
»  '  ■*•••»*:  studies.  Previously,  the  low  speed  setting  of  the  downstream  fan 
.1-  .  for  this  purpose.  The  increased  power  requirement  is  a  result  of 
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the  higher  mass  flow  in  combination  with  the  small  annulus  area  of  the  down¬ 
stream  ducting. 

A  brief  study  was  performed  to  determine  the  boundary  for  rotating 
stall  inception  on  the  isolated  rotor.  Rotating  stall  was  detected  in  two 
ways:  first  with  static  pressure  taps  on  the  outer  casing  at  the  axial  location 
of  the  rotor  quarter  chord,  and  second,  with  the  total  pressure  rake  located 
immediately  downstream  of  the  rotor.  In  the  latter  case,  the  rake  was  traversed 
over  the  full  span  of  the  rotor  blades  to  determine  the  radial  extent  of  the 
rotating  stall.  Rotating  stall  inception  was  determined  for  two  values  of  the 
mean  axial  velocity  through  the  test  section.  The  results  are  summarized  in 
Table  4. 


Table  4 

ROTATING  STALL  INCEPTION  ON  LOW 
HUB-TO-TIP  RATIO  ROTOR 


MEAN  AXIAL 
VELOCITY,  UQ,  ft/sec 

ROTOR  RPM 

AT  INCEPTION 

FLOW  COEFFICIENT 

AT  INCEPTION 

REMARKS 

33.2 

645 

0.532 

INTERMITTENT 

FULL  SPAN  STALL. 

NO  HYSTERESIS 

59.6 

1150 

0.536 

The  rotating  stall  which  occurred  on  the  low  hub-to-tip  ratio  rotor 
did  not  display  hysteresis  between  inception  with  increasing  rpm  and  dis¬ 
appearance  of  stall  with  decreasing  rpm.  Near  inception,  rotating  stall  was 
intermittent  and  was  preceded  by  gradually  increasing  turbulent  pressure  fluc¬ 
tuations  on  the  pressure  transducers  used  to  detect  rotating  stall.  This  occur¬ 
red  both  with  the  outer  casing  static  pressure  taps  and  with  the  total  pressure 
rake  in  the  annulus  downstream  of  the  rotor.  Radial  surveys  with  the  total 
pressure  rake  showed  that  rotating  stall  encompassed  the  full  span  of  the  rotor. 
Table  IV  presents  the  inception  results  measured  at  both  values  of  mean  axial 
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flow  velocity.  The  flow  coefficient,  U0/wb  ,  at  inception  for  both  mean 
axial  velocities  agree  quite  well.  This  indicates  that  inception  is  not  a 
function  of  Reynolds  number  over  the  speed  range  of  these  tests. 

The  experience  gained  in  the  preliminary  tests  on  the  modified 
annular  cascade  has  shown  that  the  rig  is  capable  of  operating  in  the  range 
required  for  studies  of  rotating  stall.  The  downstream  constriction  of  the 
annulus  used  to  accommodate  the  existing  hydraulic  drive  and  outlet  ducting, 
and  the  associated  increased  power  requirements  for  the  downstream  fan,  cause 
a  large  increase  in  operating  noise.  As  it  now  stands,  the  annular  cascade 
would  not  be  suitable  for  the  performance  of  acoustic  studies  such  as  those 
reported  in  References  13  and  14. 

Further  modifications  to  the  annular  cascade  are  planned  to  incor¬ 
porate  a  torque  measuring  apparatus  on  the  rotor  drive  shaft.  This  will 
require  disassembly  and  rework  of  the  rotor  drive  system.  In  view  of  the 
operating  noise  problem,  it  is  planned  to  remove  the  constriction  in  the 
downstream  ducting  at  the  time  the  torquemeter  is  installed.  This  should 
eliminate  the  operating  noise  problem  so  that  additional  acoustic  studies 
can  be  performed. 


E .  SUMMARY 

The  current  work  on  rotating  st^ll  in  the  annular  cascade  was 
designed  to  provide  basic  information  on  rotor-stator  interference  and  provide 
input  data  for  the  rotating  stall  stability  theory.  This  data  was  obtained 
with  the  0.8  hub-to-tip  ratio  configuration.  The  basic  experimental  data 
required  as  inputs  to  the  rotating  stall  theory  are  the  steady-state  losses 
across  the  compressor  blade  row.  These  are  required  in  a  coordinate  system 
fixed  to  the  blade  row  (relative  losses).  The  determination  of  the  relative 
losses  through  a  rotor  row  requires  the  measurement  of  total  pressure  rise 
in  a  duct -fixed  coordinate  system  and  a  measurement  of  the  work  done  by  the 
rotor.  In  previous  tests  on  an  isolated  rotor,  the  work  done  by  the  rotor  was 
determined  from  the  velocity  components  measured  downstream  of  the  rotor  with 
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a  crossed  hot-wire  anemometer  system.  In  this  case,  good  agreement  was  found 
between  the  prediction  of  the  stability  theory  and  the  experimentally  deter¬ 
mined  stall  boundary.  However,  the  situation  for  a  closely  coupled  stage  is 
more  complex  because  we  could  not  measure  the  conqponents  of  velocity  between 
the  rotor  and  stator.  Thus,  an  attempt  was  made  to  determine  the  rotor  work 
by  measuring  the  total  temperature  rise  across  the  rotor,  instead  of  the 
velocity  components.  The  results  of  the  total  temperature  studies  demonstrated 
that  the  rotor  work  coefficients  calculated  from  the  total  temperature  measure¬ 
ments  between  the  rotor  and  stator  were  reasonably  consistent  in  the  mid-annulus 
region  with  the  work  coefficients  calculated  from  velocity  measurements  down¬ 
stream  of  the  isolated  rotor.  However,  the  repeatability  of  the  total  tem¬ 
perature  measurements  was  not  good,  and  the  data  appeared  to  be  contaminated 
by  heat  addition  from  the  hydraulic  drive  for  the  rotor.  It  was  concluded 
that  the  total  temperature  results  were  not  accurate  enough  to  provide  an 
estimate  of  the  work  done  by  the  rotor  for  use  in  calculating  the  relative 
losses  across  the  rotor  in  a  rotor-stator  stage. 

On  the  other  hand,  an  extensive  series  of  total  pressure  surveys 
showed  that  at  a  given  rotor  speed,  there  was  significantly  more  total  pressure 
rise  through  the  rotor  with  the  stator  behind  it,  compared  to  the  isolated 
rotor  case.  This  could  be  due  to  either:  (a)  the  rotor-stator  interference 
produces  more  rotor  work  than  the  isolated  case,  or  (b)  the  rotor-stator  inter¬ 
ference  produces  a  reduction  in  the  rotor  losses  (presumably  by  influencing 
boundary  layer  separation) .  For  application  of  the  stability  theory,  the 
latter  assumption  was  made  and  the  work  through  the  rotor  in  the  stage  was 
assumed  to  be  the  same  as  for  the  isolated  rotor  case  at  the  same  rpm.  Thus, 
the  relative  losses  through  the  rotor  were  calculated  from  the  total  pressure 
rise  across  the  rotor  in  the  stage  and  the  work  done  by  the  isolated  rotor. 

The  rotating  stall  stability  theory  requires  several  other  inputs, 
as  well  as  the  relative  losses  through  the  rotor  in  the  rotor-stator  stage. 
These  inputs  are:  the  turning  angle  of  the  flow  through  the  rotor,  the  losses 
through  the  stator  row,  and  the  turning  angle  of  the  flow  through  the  stator 
row.  The  methods  used  to  estimate  these  parameters  from  the  experimental  data 
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are  presented  in  the  text,  along  with  the  results.  The  results  are  presented 
for  two  stagger  angle  settings  of  the  stators  in  the  stage. 

It  will  be  shown  in  the  following  section  that  the  use  of  the  above 
steady-state  data  for  input  to  the  rotating  stall  stability  theorv  leads  to 
predictions  of  rotating  stall  inception  for  the  stage  which  are  in  almost 
exact  agreement  with  the  experimental  results.  Successful  application  of  the 
theory  with  the  inputs  derived  as  discussed  above  suggests  that  rotor-stator 
interference  has  no  effect  on  the  work  done  by  the  rotor  in  the  stage. 

However,  this  is  not  conclusive  proof  as  to  which  effect  (more  rotor  work  or 
reduced  rotor  losses)  is  the  cause  of  the  significant  increase  in  total  pressure 
rise  through  the  rotor  when  the  stator  is  placed  behind  it.  Resolving  this 
point  is  an  important  factor  for  rotating  stall  research  as  well  as  for  general 
compressor  design  practice.  Hence,  it  is  planned  to  perform  additional  experi¬ 
ments  which  will  attempt  to  prove  conclusively  which  effect  is  occurring. 


The  above  experimental  studies  on  the  rotor-stator  stage  were 
performed  with  undistorted  inlet  flow.  In  addition  to  these  experiments,  a 
brief  study  was  performed  to  investigate  the  effect  of  circumferential  inlet 
distortion  on  rotating  stall  inception  in  the  rotor-stator  stage.  Here  again, 
two  stator  stagger  angles  were  used.  The  inlet  distortion  delayed  rotating 
stall  inception  by  a  small  amount  for  both  stator  stagger  angles.  Possible 
reasons  for  this  surprising  result  are  discussed  in  the  text.  Analysis  of 
the  experimental  data  to  estimate  the  inputs  required  by  the  rotating  stall 
stability  theory  is  presented. 

The  final  portion  of  this  section  describes  modifications  to  the 
annular  cascade  facility  to  reduce  the  hub-to-tip  ratio  from  0.8  to  0.44. 
Preliminary  tests  of  the  modified  rig  showed  that  it  is  capable  of  operating 
in  the  range  required  for  studies  of  rotating  stall.  The  rotating  stall 
inception  boundary  for  a  low  hub-to-tip  ratio  rotor  is  presented. 
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SECTION  IV 


COMPARISON  BETWEEN  THEORY  AND  EXPERIMENT  FOR  A 
CLOSELY  COUPLED  ROTOR-STATOR  STAGE 

A.  INTRODUCTION 

A  stability  theory  for  the  prediction  of  rotating  stall  inception  on 
a  stage  consisting  of  a  rotor  and  a  stator  was  developed  in  Reference  2.  Initial 
application  of  the  theory  indicated  that  the  addition  of  a  stator  row  down¬ 
stream  of  a  rotor  can  stabilize  the  stage,  a  result  which  was  in  agreement  with 
experimental  data  obtained  in  the  annular  cascade  facility  (Reference  1). 

However,  quantitative  comparisons  between  theory  and  experiment  for  conditions 
at  inception  of  rotating  stall  were  not  satisfactory  because  of  large  and 
questionable  extrapolations  of  the  available  steady  state  data  required  as 
inputs  to  the  theory.  The  experiments  required  to  provide  the  steady-state 
input  data  for  application  of  the  theory  are  presented  in  Section  III.  Here 
the  theoretical  predictions  will  be  compared  with  the  experimental  results. 

Comparisons  between  theory  and  experiment  for  an  isolated  rotor,  a 
closely  coupled  rotor-stator  stage,  and  the  same  stage  with  circumferential 
inlet  distortion  are  presented  in  Section  IV-B.  The  results  of  a  theoretical 
investigation  of  the  effect  of  changing  the  axial  spacing  between  the  rotor 
and  stator  in  a  stage  are  presented  in  Section  IV-C.  The  overall  results  are 
summarized  in  Section  IV-D. 

B.  COMPARISONS  BETWEEN  THEORY  AND  EXPERIMENT 

Comparisons  between  experimental  and  predicted  inception  of  rotating 

stall  on  an  isolated  rotor  were  presented  in  Reference  2  for  three  stagger  angle 

settings  of  the  rotor  (6^  =  50,  40  and  50  degrees).  The  predicted  inception 

points  agreed  very  well  with  the  experimentally  observed  inception  points  in 

all  three  cases.  However,  the  experimental  steady-state  loss  data  for  5  =  40 

KrI 

degrees  contained  a  small  error  near  rotating  stall  inception.  These  data  have 
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been  corrected  in  the  current  study.  In  addition,  the  referenced  loss  data  were 
obtained  with  a  total  pressure  rake  situated  3.9  inches  from  the  rotor  trailing 
edge.  In  the  current  program,  the  rotor  losses  were  measured  at  this  axial 
location  and  at  an  axial  location  0.7  inches  from  the  rotor  trailing  edge.  The 
results,  shown  in  Figure  14,  were  different,  with  slightly  lower  losses  calculated 
from  the  total  pressure  measurements  close  to  the  rotor.  Thus  the  stability 
theory  was  applied  using  both  loss  curves  as  inputs  to  the  theory.  Only  those 
experimental  data  measured  prior  to  rotating  stall  inception  were  used  to  estimate 
the  relative  losses  through  the  rotor.  The  lines  drawn  through  the  data  points 
in  Figure  14  are  the  inputs  used  in  the  theory.  In  addition  to  the  losses  through 
the  rotor,  the  theory  requires  the  steady-state  turning  performance  of  the  rotor. 
These  data  were  estimated  as  discussed  in  Section  III  and  are  shown  in  Figure  15. 

Computed  damping  factors  for  the  isolated  rotor  are  shown  in  Figure  21 
as  a  function  of  rotor  rpm.  The  results  using  loss  data  calculated  from 
measurements  with  the  total  pressure  rake  at  both  axial  locations  are  shown. 
Theoretical  rotating  stall  inception  occurs  when  the  damping  factor  first  goes 
to  zero.  Although  the  damping  factors  calculated  from  the  two  sets  of  loss 
measurements  differ  at  low  rotor  speeds,  they  become  the  same  at  higher  rotor 
speeds  and  indicate  the  same  rotor  speed  (998  rpm)  at  inception.  Experimentally 
observed  rotating  stall  inception,  shown  by  an  arrow  in  Figure  21,  occurred  at 
1000  rpm.  Thus  the  predicted  and  experimental  inception  points  for  the  isolated 
rotor  agree  almost  exactly  in  this  case. 

The  inputs  necessary  for  application  of  the  stability  theory  to  the 
rotor-stator  stage  require  two  other  inputs  in  addition  to  the  rotor  loss  and 
turning  performance  presented  in  Figures  14  and  15.  These  are  the  losses  through 
the  stator  row  and  the  turning  through  the  stator  row.  The  latter  inputs  were 
taken  from  Figures  16  and  17.  In  each  instance,  the  input  data  appropriate  to 
the  stator  stagger  angle  under  investigation  was  used  in  the  theory. 

Computed  damping  factors  for  the  stage  are  shown  in  Figures  22  and  23 
for  stator  stagger  angles  of  28.2  and  37.2  degrees  respectively.  In  each  of  these 
figures,  the  dashed  curve  labeled  "rotor  alone"  shows  the  damping  factor  for  the 
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rotor  without  the  stator  but  with  the  losses  measured  for  the  rotor  in  the  rotor- 
stator  stage.  The  solid  curves  in  Figures  22  and  23  show  the  damping  factors 
for  the  complete  rotor-stator  stage. 

The  first  application  of  the  theory  to  the  stage  with  a  stator  stagger 
angle  of  28.2  degrees  used  a  linear  relation  between  stator  losses  and  rotor 
speed  in  the  higher  rotor  rpm  range.  In  this  case,  the  stator  loss  curve  between 
850  and  1025  rpm  was  taken  as  the  straight  dotted  line  shown  in  Figure  16.  The 
predicted  result  for  the  stage  damping  factor  is  not  shown  in  Figure  22  because 
there  was  very  little  difference  between  this  curve  and  the  curve  for  the  rotor 
alone.  In  both  cases,  rotating  stall  inception  was  predicted  at  a  rotor  speed 
between  1015  and  1020  rpm,  whereas  experimental  inception  was  observed  at  a 
rotor  speed  of  1035  rpm.  While  the  agreement  between  these  theoretical  predictions 
and  experiment  appears  to  be  quite  good,  the  use  of  the  straight-line  stator  loss 
approximation  ignores  the  unstalled  stator  loss  measured  at  a  rotor  speed  of 
1030  rpm  (see  Figure  16).  This  last  unstalled  loss  measurement  indicates  that 
the  stator  loss  coefficient  began  to  decrease  as  rotating  stall  inception  was 
approached.  A  second  application  of  the  stability  theory  was  made  to  see  if  the 
apparent  drop  in  stator  loss  coefficient  near  stall  inception  has  a  significant 
effect  on  inception.  In  this  case,  the  straight  line  stator  loss  approximation 
was  replaced  by  a  spline  fit  to  the  unstalled  data  points  (solid  curve  in  Figure 
16).  In  cases  where  two  loss  measurements  were  made,  the  spline  fit  used  the 
average  of  the  two  data  points.  The  measurement  at  1030  rpm  was  included  in  the 
spline  fit.  With  this  stator  loss  curve,  the  calculated  stage  damping  (solid 
curve  in  Figure  22)  predicted  a  rotating  stall  inception  point  at  1042  rpm.  This 
predicted  value  is  in  slightly  better  agreement  with  the  experimental  inception 
point  (1035  rpm)  than  the  prediction  based  on  the  straight  line  loss  approxima¬ 
tion  (1017  rpm).  Thus  it  is  concluded  that  the  drop  in  stator  loss  coefficient 
just  before  rotating  stall  inception  has  a  small,  but  significant  effect  on  the 
stability  of  the  rotor-stator  stage. 

The  predicted  stability  of  the  stage  with  a  stator  stagger  angle  of 
37.2  deg.  is  shown  in  Figure  23.  In  this  case,  there  is  a  substantial  difference 
between  the  stability  of  the  rotor  alone  and  the  stability  of  the  complete  stage. 

The  rotor  alone  has  a  region  of  instability  between  rotor  speeds  of  939  and  983 
rpm  while  the  stage  is  highly  stable  in  this  region.  Rotating  stall  inception 
for  the  stage  is  predicted  at  a  rotor  speed  of  1109  rpm.  Experimental  inception 
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was  observed  at  a  rotor  speed  of  1105  rpm.  The  agreement  between  theory  and 
experiment  is  excellent  in  this  case.  It  is  worth  noting  that  for  this  stagger 
angle  (  =  37.2  deg.),  the  measured  stator  loss  curve  (Figure  16)  began  to 

drop  significantly  in  the  rotor  speed  range  between  approximately  900  and  1000 
rpm.  This  is  the  same  speed  range  where  the  rotor  alone  displays  an  unstable 
region  and  the  stage  displays  high  stability.  This  behavior  reinforces  the 
conclusion  that  the  stator  loss  characteristics  can  have  a  significant  effect 
on  the  stability  of  the  stage.  It  may  be  possible  to  take  advantage  of  the 
stabilizing  effect  of  the  stator  loss  curve  in  designing  a  compressor.  For 
instance,  it  may  be  possible  to  design  a  highly  stable  compressor  by  offsetting 
the  locations  of  the  drag  buckets  for  the  rotor  and  stator,  provided,  of  course, 
that  the  stage  efficiency  does  not  decrease  to  unacceptable  levels. 

Past  experience  with  the  two-blade-row  stability  theory  has  shown  that 
the  theory  generally  results  in  a  minimum  number  of  stall  cells  below  which  no 
unstable  solutions  are  found,  and  that  this  minimum  number  corresponds  to  the 
number  of  cells  observed  experimentally.  A  similar  result  was  found  in  this 
study;  one  stall  cell  was  found  in  the  experiments  and  one  stall  cell  was 
predicted  by  the  theoryas  the  minimum  number  of  cells  for  an  unstable  solution 
on  the  rotor-stator  stage. 

Finally,  rotating  stall  propagation  velocities  were  not  predicted  very 
well  by  the  theory.  This  is  illustrated  in  Figure  24  which  shows  the  theoretical 
propogation  velocities  relative  to  the  rotor.  The  range  of  results  from 
experimental  measurements  are  shown  as  shaded  bands  for  comparison.  The  experi¬ 
mental  measurements  showed  relative  propagation  velocites  between  35  and  45 
percent  of  the  rotor  speed  for  the  isolated  rotor.  (This  corresponds  to  absolute 
propagation  velocities  between  65  and  55  percent  of  the  rotor  speed  in  the 
direction  of  the  rotor  velocity.)  Similar  measurements  for  the  stage  lie  in  the 
range  between  approximately  60  and  70  percent  of  the  rotor  speed.  Conversely, 
the  predicted  disturbance  propagation  velocities  vary  widely  with  rotor  speed, 
particularly  for  the  stage  with  6SM  =28.2  degrees.  At  theoretical  inception, 
the  predicted  propagation  velocities  relative  to  the  rotor  are  too  large  in 
all  cases. 
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The  poor  predictions  of  stall  cell  propagation  velocity  may  result  from 
the  fact  that  the  theory  is  a  linearized  theory  and  cannot  be  expected  to  apply 
after  rotating  stall  has  started.  On  the  other  hand,  the  experimental  measure¬ 
ment  of  propagation  velocity  can  be  made  only  after  rotating  stall  has  started. 

The  theory  does  predict  one  trend  in  propagation  velocity  which  appears  con¬ 
sistent  with  the  experimental  data.  The  predicted  relative  propagation  velocities 
for  the  isolated  rotor  are  approximately  50  to  60  percent  lower  than  those  for 
the  rotor-stator  stage.  Similarly,  experimental  relative  propagation  velocities 
on  the  isolated  rotor  are  approximately  60  percent  lower  than  those  on  the  stage. 

The  final  comparison  between  theory  and  experiment  is  for  the  rotor- 
stator  stage  with  circumferential  inlet  distortion.  Experimentally,  the  dis¬ 
tortion  delayed  rotating  stall  inception  on  the  stage  for  both  stator  stagger 
angles  that  were  tested.  The  theory  was  applied  to  see  if  the  delay  in 
inception  can  be  predicted  from  the  stage  performance  measured  in  the  absence 
of  distortion  but  corrected  for  a  one  degree  decrease  in  turning  of  the  flow 
through  the  rotor.  The  justification  for  this  procedure  and  the  method  of 
data  analysis  used  to  provide  the  steady-state  inputs  to  the  theory  are  presented 
in  Section  III.C.5. 

Computed  damping  factors  for  the  stage  with  distortion  are  shown  in 
Figure  25  for  the  stage  with  a  stator  stagger  angle,  <5SM  =  28.2  degrees,  and  in 
Figure  26  for  6^37.2  degrees.  Here  again,  the  stability  of  the  complete  stage, 
and  of  the  rotor  alone  with  losses  corresponding  to  the  rotor  in  the  stage  are 
shown.  Comparison  of  these  figures  with  the  corresponding  figures  (Figures  22 
and  23)  for  undistorted  flow  shows  that  the  damping  curves  for  the  rotor  alone 
are  changed  slightly  and  those  for  the  stage  are  changed  considerably.  However, 
predicted  inception  on  the  stage  is  only  delayed  a  small  amount  in  both  cases. 

For  <S5M  =  28.2  degrees,  predicted  inception  changes  from  1042  rpm  without  dis¬ 
tortion  to  1049  rpm  with  distortion.  For  65m* 37.2  degrees,  the  corresponding 
predictions  are  1109  rpm  without  distortion  and  1113  rpm  with  distortion.  While 
the  predicted  delays  in  rotating  stall  incpetion  show  the  proper  trend,  the 
magnitudes  of  the  delays  are  not  large  enough  to  provide  accurate  predictions 
of  the  experimentally  observed  rotating  stall  inception  points.  The  latter 
occurred  at  1085  rpm  for  6Sh1  =  28.2  degrees  and  at  approximately  1160  rpm  for 
<SSm  =  37.2  degrees. 
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Application  of  the  theory  to  the  isolated  rotor  with  distortion  to 
obtain  numerical  results  has  not  been  undertaken  because  for  this  case  the 
results  are  obvious.  The  only  input  change  required  to  consider  the  effect  of 
distortion  on  the  isolated  rotor  is  a  reduction  of  the  flow  turning  through  the 
rotor  by  one  degree  and  a  corresponding  small  change  in  the  work  performed  by 
the  rotor.  The  predicted  inception  point  would  be  increased  one  or  two  rotor 
rpm's  at  most  over  that  predicted  for  undistorted  inlet  flow.  In  this  case, 
the  result  is  in  agreement  with  the  experimental  observations;  on  the  isolated 
rotor,  stationary  inlet  distortion  had  practically  no  effect  on  experimental 
rotating  stall  inception  (see  Table  3). 

The  apparent  failure  of  the  theory  to  predict  rotating  stall  inception 
accurately  for  a  stage  with  inlet  distortion  may  be  due  to  the  assumptions  made 
in  estimating  the  steady-state  blade  row  performance.  These  assumptions  were 
based  on  the  results  of  extensive  measurements  made  on  the  isolated  rotor  with 
distortion  and  may  not  be  representative  of  the  rotor  in  the  stage.  The 
performance  of  even  more  extensive  measurements  to  obtain  accurate  total  pressure 
data  for  the  rotor  in  the  stage  with  distortion  was  not  practical.  Even  with  these 
results,  it  would  still  be  necessary  to  use  the  isolated  rotor  work  for  calcu¬ 
lating  relative  losses  as  discussed  in  Section  III. 

An  alternative  possibility  for  the  inaccurate  theoretical  predictions 
of  inception  on  the  stage  with  distortion  is  the  presence  of  destructive  inter¬ 
ference  between  the  distortion  pattern  and  the  rotating  stall  which  occurs 
naturally  on  the  stage  without  distortion.  This  effect,  which  has  been  mentioned 
in  Section  III,  was  found  to  delay  inception  on  the  isolated  rotor  when  the  dis¬ 
tortion  pattern  was  rotated  about  the  compressor  axis  at  approximately  half  of 
the  rotating  stall  propagation  velocity  (Reference  1,  Figure  28).  On  the  stage, 
the  absolute  propagation  velocity  was  quite  slow,  about  30  to  40  percent  of  the 
rotor  speed.  Thus,  if  destructive  interference  occurs  on  the  stage  in  the  same 
way  as  it  did  on  the  isolated  rotor,  it  would  have  a  maximum  effect  when  the 
distortion  pattern  rotated  near  15  to  20  percent  of  the  rotor  speed.  This  is 
slow  enough  that  the  range  of  occurrence  could  encompass  stationary  distortion. 

The  magnitude  of  the  delay  in  inception  on  the  isolated  rotor  with  distortion 
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attained  a  maximum  of  about  10  percent  over  that  observed  on  the  rotor  without 
distortion.  This  is  easily  large  enough  to  explain  the  approximately  5  percent 
increase  in  rotor  speed  for  inception  on  the  stage  with  distortion. 


C.  THEORETICAL  EFFECTS  OF  BLADE  ROW  SPACING 

The  previous  comparisons  between  theory  and  experiment  demonstrated  •  t 
the  stability  theory  is  highly  accurate  in  predicting  rotating  stall  inception 
on  a  rotor-stator  stage  without  distortion.  These  comparisons  were  made  for  a 
stage  with  the  blade  rows  situated  close  together.  Some  further  theoretical 
calculations  have  been  made  to  display  the  effect  of  changing  the  axial  spacing 
between  the  rotor  and  stator  in  the  stage.  The  results  are  presented  below. 

The  steady-state  blade  row  performances  used  in  these  calculations 
were  taken  to  be  independent  of  blade  row  spacing.  That  is,  the  rotor  and 
stator  losses  and  turning  were  the  same  as  those  used  for  the  closely-coupled 
stage  discussed  previously.  Strictly  speaking,  it  is  incorrect  to  hold  the 
blade-row  performance  constant  as  the  axial  spacing  is  changed.  The  experiments 
on  the  isolated  rotor  and  on  the  rotor  in  the  stage  showed  that  losses  through 
the  rotor  (Figure  14)  change  considerably  when  the  rotor  has  a  stator  row 
added  immediately  downstream.  Thus,  one  would  expect  the  rotor  losses  to  change 
from  the  closely-coupled  stage  values  to  the  isolated  rotor  values  as  the 
distance  between  the  blade  rows  is  increased.  Nevertheless,  it  is  instructive 
to  examine  the  stability  of  the  system  as  a  direct  function  of  blade-row 
spacing,  without  considering  the  (unknown)  effects  of  spacing  changes  on  blade- 
row  performance. 

The  stage  configuration  selected  for  this  study  used  a  stator  stagger 
angle  of  37.2  degrees.  In  the  previous  results  for  this  configuration,  the 
theoretical  results  displayed  a  large  difference  between  the  stability  of  the 
rotor  alone  (but  with  the  losses  of  the  rotor  in  the  stage),  and  the  stability  of 
the  stage  (Figure  23).  The  large  differences  facilitate  comparisons  to  see  how 
the  stage  stability  approaches  that  of  the  rotor  alone  as  the  axial  spacing 

77 


/ 


between  the  blade  rows  is  increased.  The  closely-coupled  stage  used  in  the 
comparisons  between  theory  and  experiment  had  an  axial  spacing  of  1.5  inches 
between  the  mid-chords  of  the  rotor  and  stator  blades  (Figure  4).  The  additional 
theoretical  calculations  used  mid-chord  axial  spacings  of  3,  6  and  12  inches. 

The  results  are  shown  in  Figures  27(a),  (b)  and  (c) . 

Computed  damping  factors  for  the  stage  with  different  blade  row  spacings 
are  shown  in  Figures  27(a),  (b)  and  (c)  as  a  function  of  rotor  speed.  The  damp¬ 
ing  factor  of  the  rotor  alone  is  shown  for  comparison.  Each  part  of  Figure  27 
shows  the  calculated  damping  for  a  given  number  of  stall  cells  and  variable 
axial  spacing.  Figure  27(a)  presents  the  results  for  one  stall  cell,  27(b) 
for  three  stall  cells,  and  27(c)  for  seven  stall  cells.  Larger  numbers  of 
stall  cells  provided  results  very  close  to  those  found  for  seven  cells  (Figure 
27(c)).  For  the  rotor  alone,  the  magnitude  of  the  damping  factor  is  a  function 
of  the  number  of  stall  cells,  but  the  rotor  speeds  where  the  damping  equals 
zero  are  not  a  function  of  number  of  cells.  Thus,  the  regions  of  instability 
for  the  rotor  alone  remain  unchanged  in  all  three  parts  of  Figure  27.  However, 
for  the  stage,  the  zero  damping  points  do  depend  on  the  number  of  stall  cells. 
This  can  be  seen  by  comparing  the  stage  damping  factors  in  Figures  27(a),  (b) 
and  (c) . 

In  the  experiments  on  the  close-coupled  stage  (T  =  1.5  inches), 
rotating  stall  inception  was  found  to  occur  at  a  rotor  speed  of  1105  rpm  with 
one  stall  cell  propagating.  It  was  noted  in  the  comparison  between  theory  and 
experiment  for  this  case  that  the  theory  generally  results  in  a  minimum  number 
of  stall  cells  below  which  no  unstable  solutions  are  found  and  that  this  minimum 
number  corresponds  to  the  number  of  cells  observed  experimentally.  This  result 
can  be  seen  by  comparing  the  damping  factors  for  Tc  =  1.5  inches  in  Figure  27(a), 
(b)  and  (c).  With  one  stall  cell  (Figure  27(a)),  the  damping  factor  for  this 
case  becomes  zero  at  1109  rpm.  The  corresponding  rotor  speeds  for  higher  numbers 
of  cells  are  1112  rpm  with  three  cells  (Figure  27(b))  and  1127  rpm  with  seven 
cells  (Figure  27(c)).  Thus  theoretical  rotating  stall  inception  occurs  first 
at  1109  rpm  with  one  stall  cell.  This  result  is  in  excellent  agreement  with  the 
experiments. 
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Now  consider  what  happens  as  the  blade  row  spacing  is  increased.  For 

one  stall  cell  (Figure  27(a)),  increasing  the  spacing  has  little  effect  on 

predicted  inception  (1109  to  1111  rpm) .  For  three  stall  cells  (Figure  27(b)), 

increased  spacing  makes  the  stage  slightly  more  stable  at  first  (T\  =  5  inchesj 

and  then  progressively  less  stable,  approaching  that  of  the  rotor  alone  at 

=  12  inches.  Kith  seven  stall  cells  (Figure  27(c)),  the  damping  factor  for 

T  =  6  inches  is  close  to  that  of  the  rotor  alone  and  for  T  =12  inches  it  is 
c  c 

the  same  as  that  of  the  rotor  alone.  Overall  inspection  of  the  results  in 
Figures  2"  (a),  (b)  and  (c)  shows  that  the  stage  stability  is  a  function  of  both 
the  axial  spacing  between  the  blade  rows  and  the  number  of  rotating  stall 
cells.  Table  5  summarises  the  unstable  rotor  speed  ranges  as  a  function  of 
blade  row  spacing  and  number  of  stall  cells.  For  the  rotor  alone  and  for  6  and 
12  inch  blade-row  spacings  in  the  stage,  there  is  an  instability  at  intermediate 
rotor  speeds  as  well  as  a  final  instability  at  a  higher  rotor  speed.  Both 
unstable  ranges  are  given  in  the  table.  The  numbers  of  stall  cells  shown  in 
each  range  were  taken  from  the  criterion  that  the  minimum  number  predicted 
theoretically  will  usually  occur.  It  is  evident  that  the  stability  of  the 
system  changes  with  blade-row  spacing  and  that  the  number  of  rotating  stall 
cells  is  a  nonlinear  function  of  the  spacing. 

The  above  theoretical  results  display  the  effect  of  axial  spacing  in 
a  stage  for  the  situation  where  the  spacing  does  not  affect  the  steady-state 
performance  of  the  blade  rows.  As  discussed  earlier,  there  is  experimental 
evidence  that  the  steady-state  blade-row  performance  is  not  independent  of  the 
axial  separation  between  the  blade  rows.  Thus  the  assumption  of  constant  blade- 
row  performance  is  not  correct,  and  the  theoretical  predictions  should  be  used 
only  as  a  guide  to  indicate  trends.  In  this  respect,  the  calculations  show 
that  changes  in  blade-row  spacing  change  both  the  stability  of  the  stage  and 
the  number  of  stall  cells  which  occur  at  inception.  The  theory  suggests  that 
larger  blade-row  spacings  will  decrease  the  stability  of  the  stage  and  increase 
the  number  of  stall  cells  which  occur  at  inception.  However,  both  the  stability 
and  number  of  cells  are  nonlinear  functions  of  blade-row  spacing. 
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1095  SEVEN  AT  INCEPTION,  PROBABLY 

CHANGING  TO  THREE  THEN  ONE. 


D. 


SUMMARY 


* 


The  experimental  steady-state  measurements  on  an  isolated  rotor  and  on 
a  rotor-stator  stage  have  been  used  as  inputs  to  the  rotating  stall  stability 
theory.  Final  application  of  the  theory  to  the  stage  configuration  required  the 
use  of  data  measured  on  the  isolated  rotor  as  well  as  data  measured  on  the  stage. 
The  theoretical  predictions  of  rotating  stall  inception  and  of  the  number  of 
stall  cells  which  occur  were  found  to  agree  extremely  well  with  the  experimental 
results.  However,  prediction  of  stall  cell  propagation  velocity  was  poor.  The 
theoretical  results  suggest  that  stator  losses  as  well  as  the  rotor  losses  play 
a  major  role  in  determining  the  stage  stability.  It  may  be  possible  to  design 
a  compressor  which  takes  advantage  of  the  stator  loss  characteristics  to  increase 
compressor  stability. 

A  further  application  of  the  theory  was  made  to  see  if  the  effects  of 
inlet  distortion  could  be  predicted.  Input  data  for  the  theory  were  estimated 
from  steady-state  data  obtained  on  the  isolated  rotor  with  distortion  and  from 
stage  data  without  distortion.  For  the  isolated  rotor,  both  the  theory  and  the 
experiments  showed  negligible  effects  of  inlet  distortion.  For  the  stage,  the 
theory  predicted  rotating  stall  inception  would  be  delayed  only  a  very  small 
amount  by  inlet  distortion,  while  the  experiments  displayed  a  larger  delay  in 
inception. 

A  final  application  of  the  theory  was  made  to  study  the  theoretical 
effect  of  changing  the  axial  separation  between  the  blade  rows  in  a  stage.  For 
this  theoretical  study,  it  was  assumed  that  blade-row  axial  spacing  has  no 
effect  on  the  steady-state  blade-row  performance.  The  experimental  data  showed 
that  this  assumption  is  not  correct.  Nevertheless,  its  use  allowed  a  study  of 
the  direct  effect  of  geometry  on  the  stage  stability.  The  results  suggest  that 
large  blade-row  spacings  will  decrease  the  stability  of  the  stage  and  increase 
the  number  of  stall  cells  which  occur  at  inception.  However,  both  the  stage 
stability  and  number  of  cells  are  nonlinear  functions  of  blade-row  spacing. 


SECTION  V 

SUMMARY  AND  CONCLUSIONS 


A  combined  experimental  and  theoretical  research  program  on  rotating 
stall  in  axial  flow  compressors  has  been  conducted.  The  theoretical  studies 
were  concerned  with  the  development  of  an  incompressible  three-dimensional 
rotating  stall  stability  theory  and  with  revision  of  a  two-dimensional  stability 
theory  for  compressible  flow.  The  major  portion  of  the  experimental  program 
was  oriented  towards  studying  rotor-stator  interference  in  a  closely  coupled 
rotor-stator  stage.  The  experimental  results  have  been  compared  with  the 
predictions  of  an  existing  two-dimensional  rotating  stall  stability  theory  for 
a  rotor-stator  stage. 

In  the  theoretical  studies,  the  two-dimensional  stability  theory  for 
wholly  subsonic  compressible  flows  was  modified  through  the  revision  of  the 
energy  matching  conditions.  Also  the  formulation  of  the  theory  to  supersonic 
flow  has  been  initiated  and  some  analysis  of  the  theory  has  been  made  with 
respect  to  the  mechanisms  of  rotating  stall  in  supersonic  flows.  In  addition 
to  the  above,  a  stability  theory  for  three-dimensional,  incompressible  flow  was 
developed  but  not  implemented  because  of  extensive  computational  requirements. 

It  appears  that  the  stability  of  three-dimensional  flows  would  be  best  studied 
through  analysis  of  the  response  of  the  blade  row  rather  than  a  stability  theory 
approach. 


The  experimental  studies  on  rotor-stator  interference  were  performed 
in  the  Calspan/Air  Force  annular  cascade  in  a  configuration  with  a  hub-to-tip 
ratio  of  0.8.  The  investigation  included  measurements  on  an  isolated  rotor,  a 
closely  coupled  rotor-stator  stage,  and  the  same  stage  with  circumferential 
inlet  distortion.  The  experiments  on  the  stage  were  performed  for  two  different 
stagger  angle  settings  of  the  stators.  With  undistorted  flow,  it  was  found 
that  the  presence  of  the  stators  delayed  rotating  stall  inception  over  that 
which  occurs  on  the  isolated  rotor.  Moreoever,  prior  to  rotating  stall  inception, 
there  was  significantly  more  total  pressure  rise  across  the  rotor  with  the  stator 
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behind  it  than  there  was  across  the  isolated  rotor.  This  could  be  due  to  either 
(a)  the  rotor-stator  interference  produces  more  rotor  work  than  the  isolated 
case,  or  (b)  the  rotor- stator  interference  produces  a  reduction  in  the  rotor 
losses  (presumably  by  influencing  boundary  layer  separation) .  An  attempt  was 
made  to  resolve  these  possibilities  by  measuring  the  total  temperature  rise 
across  the  rotor.  However,  the  results  were  not  accurate  enough  to  provide 
reliable  estimates  of  the  work  done  by  the  rotor. 

The  experimental  steady-state  performance  data  measured  on  the  isolated 
rotor  and  on  the  stage  were  used  to  provide  inputs  to  a  two-dimensional  rotating 
stall  stability  theory.  The  inputs  required  by  the  theory  are  the  relative 
losses  and  the  turning  of  the  flow  through  the  rotor  and  the  losses  and  turning 
of  the  flow  through  the  stators.  Calculation  of  the  relative  losses  through  the 
rotor  in  the  stage  requires  determination  of  the  total  pressure  rise  across  the 
rotor  and  of  the  work  done  by  the  rotor.  As  noted  above,  attempts  to  measure 
the  rotor  work  through  total  temperature  measurements  were  not  successful. 
However,  the  work  done  by  the  isolated  rotor  had  been  determined  previously 
through  an  extensive  series  of  velocity  measurements.  Thus  for  application  of 
the  stability  theory,  it  was  assumed  that  the  work  done  by  the  rotor  in  the 
stage  was  the  same  as  that  done  by  the  isolated  rotor.  This  corresponds  to 
assuming  that  the  observed  increase  in  total  pressure  rise  across  the  rotor 
in  the  stage  is  a  result  of  reduced  losses  through  the  rotor.  With  this  input, 
the  resulting  stability  boundary  predictions  from  the  theory  agreed  extremely 
well  with  the  experimental  stall  inception  boundary  for  the  stage  for  both 
stator  stagger  angles  that  were  tested.  However,  this  result  is  not  conclusive 
proof  as  to  which  effect,  more  rotor  work  or  reduced  rotor  losses,  is  the  cause 
of  the  significant  increase  in  total  pressure  rise  across  the  rotor  when  the 
stator  is  placed  behind  it.  Resolving  this  point  is  an  important  factor  for 
rotating  stall  research  as  well  as  for  general  compressor  design  practice. 

In  addition  to  predicting  the  conditions  for  rotating  stall  inception, 
the  theory  predicts  the  number  of  stall  cells  which  will  occur  and  the  propaga¬ 
tion  velocity  of  the  stall  cells.  The  predicted  number  of  stall  cells  agreed 
with  the  experiments  but  the  prediction  of  stall  cell  propagation  velocity  was 
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poor.  The  poor  predictions  of  propagation  velocity  may  result  from  the  fact  that 
the  theory  is  is  a  linearized  theory  and  cannot  be  expected  to  apply  after  rotating 
stall  has  started.  On  the  other  hand,  experimental  measurements  of  propagation 
velocity  can  be  made  only  after  rotating  stall  has  started.  The  theoretical 
results  suggest  that  the  stator  loss  characteristics  can  have  a  significant 
effect  on  the  stability  of  the  stage.  It  may  be  possible  to  take  advantage  of 
the  stabilizing  effect  of  the  stator  loss  curve  in  designing  a  compressor. 

Finally,  the  theoretical  effect  of  changing  the  axial  separation  between  the 
rotor  and  stator  in  a  stage  is  presented. 

The  above  experimental  studies  were  performed  with  undistorted  inlet 
flow.  In  addition,  a  brief  study  was  performed  to  investigate  the  effect  of 
circumferential  inlet  distortion  on  rotating  stall  inception  in  the  rotor- 
stator  stage  with  two  different  settings  of  the  stator  stagger  angle.  The 
results  were  surprising;  the  presence  of  inlet  distortion  caused  a  delay  in 
rotating  stall  inception  for  both  stator  stagger  angles.  A  further  application 
of  the  stability  theory  was  made  to  see  if  these  results  could  be  predicted. 

Input  data  for  the  theory  were  estimated  from  steady-state  data  measured  on  the 
isolated  rotor  with  distortion  and  from  stage  data  without  distortion.  With 
these  inputs  the  theory  predicted  a  delay  in  inception  but  the  delay  was  not  as 
large  as  that  observed  experimentally.  In  contrast,  for  the  isolated  rotor, 
both  the  theory  and  experiments  showed  negligible  effects  of  inlet  distortion. 

The  possibility  that  destructive  interference  between  rotating  stall  and  inlet 
distortion  caused  the  delayed  inception  on  the  stage  is  discussed  in  the  text. 

Such  a  phenomenon  has  been  observed  previously  on  an  isolated  rotor  with  a  dis¬ 
tortion  pattern  rotating  about  the  compressor  axis. 

At  the  conclusion  to  the  above  studies,  the  annular  cascade  was 
modified  to  reduce  the  hub-to-tip  ratio  of  the  annulus  from  0.8  to  0.44.  These 
modifications  are  described  along  with  the  results  of  checkout  tests  and  a 
preliminary  study  to  investigate  the  rotating  stall  boundary  on  a  low  hub-to-tip 
ratio  rotor. 
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APPENDIX  A 

FUNCTIONS  APPEARING  IN  THE  THREE-DIMENSIONAL  STABILITY  THEORY 


Eq.  (53)  appearing  in  Section  II-C  is  given  below  along  with  the 
expression  for  the  coefficients  of  the  derivatives. 
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Figure  2  OVERALL  VIEW  OF  ANNULAR  CASCADE  FACILITY 
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Figure  3  VIEW  OF  ANNULAR  CASCADE  TEST  SECTION 
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Figure  5  NOTATION  FOR  ANNULAR  CASCADE  WITH  ROTOR-STATOR  STAGE 
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OVERALL  TOTAL  PRESSURE  RISE  COEFFICIENT.  ACpj 


Figure  7  OVERALL  TOTAL  PRESSURE  RISE  ACROSS  ISOLATED  ROTOR  AND  ROTOR 
IN  ROTOR-STATOR  STAGE 
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Figure  8  RADIAL  DISTRIBUTIONS  OF  TOTAL  PRESSURE  COEFFICIENT  DOWNSTREAM 
OF  ISOLATED  ROTOR.  TOTAL  PRESSURE  RAKE  LOCATED  AT  AXIAL 
DISTANCE  3.9  INCHES  FROM  ROTOR  TRAILING  EDGE 


Figure  9  RADIAL  DISTRIBUTIONS  OF  TOTAL  PRESSURE  COEFFICIENT  DOWNSTREAM 
OF  ISOLATED  ROTOR.  COMPARISON  OF  MEASUREMENTS  CLOSE  TO  ROTOR 
WITH  MEASUREMENTS  FARTHER  DOWNSTREAM 
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Figure  11  TOTAL  PRESSURE  RISE  ACROSS  ISOLATED  ROTOR  AND  ACROSS  ROTOR 
IN  CLOSELY  COUPLED  ROTOR-STATOR  STAGE 
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ROTOR  STAGGER  ANGLE.  6BM  “  40  de 


Figure  13  RADIAL  DISTRIBUTIONS  OF  WORK  COEFFICIENT  MEASURED  ON  ISOLATED  ROTOR 


OVERALL  TOTAL  PRESSURE  LOSS  COEFFICIENT,  f-A  C„  1  = - AC 


Figure  14  OVERALL  TOTAL  PRESSURE  LOSS  THROUGH  ROTOR  IN  ROTOR-STATOR 
STAGE,  COORDINATE  SYSTEM  RELATIVE  TO  ROTOR 
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ABSOLUTE  AND  RELATIVE  OVERALL  OUTLET  SWIRLS  FROM  ISOLATED  ROTOR 


FROM  ROTOR,  tan  0 


□  MEASURED  DIRECTLY (UNSTALLED) 
cf  MEASURED  DIRECTLY  (STALLED) 
-  CALCULATED  FROM  ROTOR  WORK 
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Figure  15  COMPARISON  OF  OVERALL  OUTLET  SWIRL  FROM  ISOLATED  ROTOR 
CALCULATED  FROM  ROTOR  WORK  WITH  DIRECT  MEASUREMENTS 
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Figure  16  OVERALL  TOTAL  PRESSURE  LOSS  THROUGH  STATORS  IN  ROTOR-STATOR  STAGE 
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Figure  18  EFFECT  OF  DISTORTION  SCREEN  ON  TOTAL  PRESSURE  RISE  OF  ISOLATED 
ROTOR,  ABSOLUTE  COORDINATE  SYSTEM 
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Figure  19  EFFECT  OF  DISTORTION  SCREEN  ON  AVERAGE  SWIRL  ANGLES  DOWNSTREAM 
OF  ISOLATED  ROTOR,  ABSOLUTE  COORDINATE  SYSTEM 
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Figure  21  THEORETICAL  STABILITY  CHARACTERISTICS  OF  ISOLATED  ROTOR 


Figure  22  THEORETICAL  STABILITY  CHARACTER ISTI 
STATOR  STAGGER  ANGLE,  5SM  =  28.2  deg 
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Figure  23  THEORETICAL  STABILITY  CHARACTERISTICS  OF  ROTOR-STATOR  STAGE, 
STATOR  STAGGER  ANGLE,  5  oM  =  37.2  DEG 


THEORETICAL 

INCEPTION 


Figure  24  THEORETICAL  ROTATING  STALL  PROPAGATION  VELOCITIES 


Figure  25  THEORETICAL  STABILITY  CHARACTERISTICS  OF  ROTOR-STATOR  STAGE 
WITH  CIRCUMFERENTIAL  INLET  DISTORTION,  STATOR  STAGGER  ANGLE, 
^  cm  =  28.2  DEG 


THEORETICAL  STABILITY  CHARACTERISTICS  OF  ROTOR-STATOR  STAGE  WITH 
CIRCUMFERENTIAL  INLET  DISTORTION,  STATOR  STAGGER  ANGLE,  6  =  37.3 


AXIAL  DISTANCE  BETWEEN 
ROTOR  AND  STATOR  MIDCHORDS 

T_  (in.) 


Figure  27  EFFECT  OF  BLADE  ROW  SPACING  ON  THEORETICAL  DAMPING  FACTORS  FOR 
ROTOR-STATOR  STAGE,  5RM  =  40OEG,  6SM  =  37.2DEG 

a.  NUMBER  OF  CELLS,  n  =  1 


AXIAL  DISTANCE  BETWEEN 
ROTOR  AND  STATOR  MIDCHORDS 


27  (Cont.)  EFFECT  OF  BLADE  ROW  SPACING  ON  THEORETICAL  DAMPING  FACTORS 
FOR  ROTOR-STATOR  STAGE,  6  dM  =  40  DEG,  5  cm  =  37.2  DEG 


